Environmental  Laboratory  erdc/el CR-13-1 


US  Army  Corps 
of  Engineers® 

Engineer  Research  and 
Development  Center 


INNOVATIVE  SOLUTIONS 
for  a  safer,  better  world 


Aquatic  Plant  Control  Research  Program 

Golden  Algae  Control:  Efficacy  of  Hydraulic 
Manipulations  in  Coves  of  Lake  Granbury 

Daniel  L.  Roelke,  Bryan  W.  Brooks,  James  P.  Grover,  April  2013 

and  Veronica  M.  Lundgren 


Approved  for  public  release;  distribution  is  unlimited. 


The  US  Army  Engineer  Research  and  Development  Center  (ERDC)  solves  the 
nation’s  toughest  engineering  and  environmental  challenges.  ERDC  develops  innovative 
solutions  in  civil  and  military  engineering,  geospatial  sciences,  water  resources,  and 
environmental  sciences  for  the  Army,  the  Department  of  Defense,  civilian  agencies,  and 
our  nation’s  public  good.  Find  out  more  at  www.erdc.usace.armv.mil. 

To  search  for  other  technical  reports  published  by  ERDC,  visit  the  ERDC  online  library 
at  http://acwc.sdp.sirsi.net/client/default. 


Aquatic  Plant  Control  Research  Program 


ERDC/EL  CR-13-1 
April  2013 


Golden  Algae  Control:  Efficacy  of  Hydraulic 
Manipulations  in  Coves  of  Lake  Granbury 


Daniel  L.  Roelke,  Veronica  M.  Lundgren 

Texas  A&M  University,  Texas  AgriLife  Research 
Department  of  Wildlife  and  Fisheries  Sciences 
2258  TAMUS 
College  Station,  TX  77843 

Bryan  W.  Brooks 

Baylor  University 

Center  for  Reservoir  and  Aquatic  Systems  Research 
One  Bear  Place  ff97266 
Waco,  TX  76798-7266 

James  P.  Grover 

University  of  Texas  at  Arlington 
Department  of  Biology 
Box  19498 

Arlington,  TX  76019-0001 


Final  report 

Approved  for  public  release;  distribution  is  unlimited. 


Prepared  for  U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 

Under  Work  Unit  33143 

Monitored  by  U.S.  Army  Engineer  Research  and  Development  Center 
3909  Flails  Ferry  Road,  Vicksburg,  MS  39180-6199 


ER  DC/EL  CR-13-1 


ii 


Abstract 

Bloom  mitigation  approaches  employing  deepwater  flushing  show  great 
promise  when  targeting  coves  of  Lake  Granbury,  USA.  They  may  be  limited, 
however,  based  on  the  natural  water  exchange  rate  between  targeted  coves 
and  the  surface  waters  of  the  open  lake.  This  potential  limitation  was 
investigated  by  conducting  dye-tracer  studies  that  enabled  direct  measure¬ 
ment  of  water  exchange  rates.  Historical  data  were  also  analyzed  to  explore 
the  dissimilarity  between  coves  and  the  open  lake.  Finally,  a  numerical 
model  was  built  to  estimate  P.  parvum  bloom  reduction  levels  under 
various  deepwater  flushing  scenarios.  Dye  tracer  studies  involved  time- 
series,  high-resolution  spatial  mapping  of  Rhodamine  WT  in  three  target 
coves.  To  complement  spatial  mapping,  vertical  profiling  was  performed  at 
fixed-location  stations  throughout  these  coves.  These  combined  processes 
allowed  researchers  to  calculate  the  total  amount  of  Rhodamine  WT  over 
time  in  these  coves.  From  this  information,  natural  water  exchange  rates 
between  coves  and  surface  waters  of  the  open  lake  could  be  determined,  and 
were  typically  found  to  be  in  the  range  of  ~o. 25-0.30  d1.  Analysis  of 
historical  data  utilized  an  extensive  data  set  comprised  of  biotic  and  abiotic 
parameters  from  10  coves  and  10  open  lake  stations.  The  data  record 
spanned  multiple  years  with  a  once-per-month  sampling  frequency. 
Dissimilarity  relationships  (estimated  from  Euclidean  distance  measures) 
between  coves  and  the  open  lake  with  changes  in  winds  and  inflows  to  the 
lake  were  examined.  It  was  found  that  water  exchanges  with  coves  located 
lower  in  the  lake  were  linked  to  winds  and  inflows.  Mechanisms  driving 
water  exchanges  with  coves  higher  in  the  lake,  however,  are  still  not  clear. 
Proposed  management  scenarios  were  depicted  with  deepwaters  being 
pumped  from  the  open  lake  to  the  upper  regions  of  coves.  Results  suggest 
that  P.  parvum  blooms  can  effectively  be  mitigated  in  cove  habitats  of  Lake 
Granbury  through  managed  deepwater  hydraulic  flushing.  In  some 
scenarios,  blooms  were  reduced  by  >70%.  The  next  steps  in  research  aimed 
at  further  developing  this  mitigation  approach  should  be  a  coupled  study 
involving  model  refinement  and  a  valuation  study. 
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Preface 

This  technical  report  was  developed  to  convey  the  results  of  in-lake  tracer 
experiments,  historical  data  analyses,  and  numerical  modeling  conducted 
as  part  of  a  project  titled  “Golden  algae  control:  Efficacy  of  hydraulic 
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1  Introduction 

Background 

Harmful  algal  blooms  (HABs)  have  increased  in  frequency,  duration,  and 
magnitude  worldwide  (Smayda  1990,  Hallegraeff  1993)  and  produce 
deleterious  effects  ranging  from  fish  kills  to  human  health  risks  (Van  Dolah 
et  al.  2001).  Unfortunately,  environmental  conditions  leading  to  HABs  are 
complex  and  often  species-specific,  making  it  difficult  to  develop  a  universal 
approach  to  management  (Anderson  and  Garrison  1997,  Roelke  2000, 
Roelke  and  Buyukates  2001,  Maso  and  Garces  2006). 

The  haptophyte  Prymnesium  parvum  (golden  algae)  is  an  HAB-causing 
species,  which  occurs  worldwide  and  is  responsible  for  large  fish  kills  in 
coastal  and  inland  water  bodies  (Moestrup  1994,  Edvardsen  and  Paasche 
1998).  Recent  studies  suggest  mechanisms  by  which  P.  parvum  might  gain 
a  selective  advantage  over  other  phytoplankton,  thereby  allowing  bloom 
formation.  One  such  mechanism  is  the  production  of  toxins  released  into 
the  water  column  where  they  impact  fish,  phytoplankton,  zooplankton,  and 
bacteria  (Ulitzur  and  Shilo  1964,  Igarashi  et  al.  1995,  Nygaard  and  Tobiesen 
1993,  Tillmann  2003,  Fistarol  et  al.  2003,  Rosetta  and  McManus  2003, 
Barreiro  et  al.  2005,  Baker  et  al.  2007,  2009).  Some  of  the  toxins  apparently 
act  as  allelopathic  substances  that  inhibit  the  growth  of  competing 
phytoplankton,  and  also  act  as  grazer  inhibitors  (Fistarol  et  al.  2003; 
Sopanen  et  al.  2006;  Roelke  et  al.  2007,  2010;  Brooks  et  al.  2010). 

A  P.  parvum  bloom  was  confirmed  along  the  Pecos  River,  Texas  (USA)  in 
1985  (James  and  De  La  Cruz  1989).  Subsequently,  P.  parvum  blooms  in 
southcentral  USA  have  affected  over  21  lakes  in  five  river  basins,  killing 
millions  of  sport  fish  and  resulting  in  millions  of  dollars  in  economic 
impact.  In  Texas,  the  majority  of  the  P.  parvum  damage  occurred  over  the 
last  10  years  (Southard  et  al.  2010,  Roelke  et  al.  2011).  Nationwide,  P. 
parvum  is  now  documented  in  17  states,  mostly  in  the  southern  regions 
(Sager  et  al.  2008,  Roelke  et  al.  2011),  but  with  toxic  blooms  also  now 
occurring  in  northern  regions  (Brooks  et  al.  2011). 

Although  factors  controlling  P.  parvum  blooms  in  lakes  of  the  southcentral 
United  States  are  complex  and  partially  unknown,  the  research  documented 
in  this  report  has  significantly  advanced  understanding  of  P.  parvum 


ER  DC/EL  CR-13-1 


2 


dynamics  in  inland  waters.  However,  there  appears  to  be  a  link  between 
eutrophication,  salinization,  and  the  appearance  of  P.  parvum  blooms. 
Blooms  in  Europe  (Holdway  et  al.  1978,  Kaartvedt  et  al.  1991,  Amsinck  et  al. 
2005),  the  Middle  East  (van  Rijn  and  Shilo  1989),  and  Asia  (Krasnotshchek 
and  Abramowitsch  1971,  Guo  et  al.  1996)  have  all  occurred  in  aquatic 
systems  that  were  eutrophic  and  brackish.  P.  parvum  blooms  in  the 
southcentral  United  States  have  appeared  only  in  lakes  with  a  higher  salt 
content  (>  1  psu)  (Texas  Parks  and  Wildlife  Department  (TPWD)  2003). 
Although  P.  parvum  is  present  in  lakes  with  lower  salt  content,  no  blooms 
have  been  recorded. 

Over  the  past  few  years,  the  research  team  has  closely  examined  P.  parvum 
dynamics  in  Lake  Granbury,  an  impoundment  on  the  Brazos  River  in 
central  Texas.  The  Lake  Granbury  fishery  has  experienced  significant 
impacts  by  P.  parvum  blooms,  resulting  in  a  societal  impetus  for  develop¬ 
ment  of  viable  environmental  management  approaches  in  a  timely  fashion. 
Fortunately,  relatively  recent  observations  have  provided  information  that 
may  allow  for  the  development  of  P.  parvum  mitigation  options  in  Lake 
Granbury  and  other  inland  waters.  These  include  hydraulic  flushing  with 
lake  deepwaters  (Hayden  et  al.  2012,  Roelke  et  al.  2012a). 

The  research  team  identified  that  hydraulic  flushing  and  its  associated 
nutrient  loading  in  Lake  Granbury  terminate  system-wide,  fish-killing  P. 
parvum  blooms  (Roelke  et  al.  2010,  2011).  The  research  team  also  identified 
that  hydraulic  dilution  (in  the  absence  of  flushing)  and  its  associated 
nutrient  loading  in  Lake  Whitney  terminated  a  fish-killing  P.  parvum 
bloom  in  April  2007  (Schwierzke-Wade  et  al.  2011).  Subsequently,  it  was 
hypothesized  that  targeted  management  of  reservoir  inflows  and  hydraulic 
flushing  rates  could  mitigate  P.  parvum  blooms. 

A  limitation  to  this  approach,  however,  would  be  available  water  upstream 
of  the  impacted  area.  To  avoid  this,  it  was  suggested  that  deepwaters  of  a 
lake,  free  of  P.  parvum  cells,  could  be  pumped  to  the  surface  to  affect 
blooms  in  the  same  way  that  external  inflows  do.  This  notion  was  success¬ 
fully  tested  using  in-lake  mesocosm  experiments  (Hayden  et  al.  2012, 

Roelke  et  al.  2012a),  a  project  funded  by  the  US  Army  Corps  of  Engineers 
(USACE)  in  2010-11.  Briefly,  it  was  found  that  P.  parvum  cell  density  was 
reduced  by  69%  and  ambient  toxicity  was  completely  ameliorated  during 
the  period  of  bloom  initiation  in  the  lake.  During  conditions  of  bloom 
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maturation,  population  density  was  reduced  by  53%,  toxicity  by  57%,  and 
bloom  proportions  were  never  reached. 

This  deep-water  pumping  had  minimal  effect  on  phytoplankton  and 
zooplankton,  and  little  effect  on  water  quality.  Shifts  towards  more  rapidly 
growing  phytoplankton  taxa  were  observed,  as  was  enhanced  abundance 
of  copepod  nauplii.  In  other  words,  while  inflows  using  deep  waters 
suppressed  P.  parvum  bloom  initiation  and  development,  they  were 
benign  to  other  aspects  of  the  lower  food  web  and  environment,  an 
appealing  outcome  for  a  potential  management  approach  (Hayden  et  al. 
2012,  Roelke  et  al.  2012a). 

While  bloom  mitigation  approaches  employing  deepwater  flushing  show 
great  promise  when  targeting  coves  of  a  lake,  they  may  be  limited  based  on 
the  natural  surface  water  exchange  rate  between  targeted  coves  and  the 
mainstem  of  a  lake.  For  example,  if  the  natural  surface  water  exchange  rate 
is  too  great,  the  effects  of  deepwater  flushing  of  a  cove  may  be  masked.  If 
natural  surface  water  exchange  rates  are  not  excessive,  however,  then 
hydraulic  manipulations  in  selected  coves  might  lead  to  creation  of  refugia 
from  P.  parvum  blooms.  This  would  greatly  accelerate  recovery  of  flora  and 
fauna  in  the  lake  after  P.  parvum  blooms  subside,  including  economically 
important  fisheries  in  the  region. 

Objective  and  approach 

The  narrative  below  describes  results  from  a  project  that  addressed  the 
issue  of  water  exchange  between  coves  and  the  mainstem  of  a  lake.The 
objective  of  this  project  was  to  identify  a  cove  selection  process  in  which  in¬ 
lake  demonstrations  of  P.  parvum  bloom  mitigation  could  be  conducted. 

Project  components 

The  project  was  accomplished  by: 

1.  Directly  measuring  water  exchange  rates  between  coves  and  natural 
surface  waters  of  the  open  lake,  and  exploring  relationships  between  water 
exchange,  wind  speed  and  direction,  and  inflow  to  the  lake. 

2.  Determining  seasonal  differences  between  plankton  community 
composition  and  water  quality  observed  in  coves  and  open  lake  waters 
based  on  wind  speed  and  direction,  and  inflow  to  the  lake. 
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3.  Determining  the  sensitivity  of  P.  parvum  bloom  formation  under 

conditions  of  managed  deepwater  flushing  of  coves  to  the  rate  of  natural 
surface  water  exchange  between  coves  and  open  waters  of  the  lake. 

Relevance 

P.  parvum  blooms  are  rapidly  spreading,  and  their  proliferation  is  antici¬ 
pated  to  get  worse  with  climate  change  and  human  population  growth 
(Roelke  et  al.  2011,  2012b).  This  research  furthered  understanding  of  a 
potential  P.  parvum  mitigation  strategy  that  would  involve  hydraulic 
manipulations  of  coves,  e.g.  deepwater  hydraulic  flushing.  As  mentioned 
above,  this  approach  was  successfully  tested  in  enclosure  experiments 
conducted  within  the  lake  (Hayden  et  al.  2012,  Roelke  et  al.  2012a).  An 
issue  of  uncertainty  facing  further  development  of  this  mitigation  strategy 
must  address  uncertainty  regarding  the  duration  that  cove  treatments 
would  remain  effective.  This  duration  is  linked  to  the  rate  at  which  cove 
waters  naturally  exchange  with  surface  waters  of  the  open  lake.  With  this 
knowledge,  coves  targeted  for  creation  of  refuge  habitat  can  better  be 
selected. 

Methods 

Study  region 

Lake  Granbury  is  a  reservoir  on  the  Brazos  River,  Texas,  USA,  constructed 
in  1969.  The  lake  has  a  capacity  of  167  x  106  m3,  a  surface  area  of  34  km2, 
and  an  average  depth  of  -5  m.  The  shoreline  follows  the  meandering  river 
channel  with  an  elongated,  sinuous  basin  oriented  northwest  to  southeast 
that  is  -45  km  long  and  has  an  average  width  of  0.6  km.  Daily  discharges 
from  the  Brazos  River  into  the  lake  are  measured  at  a  location  upstream 
from  the  lake  (Dennis,  Texas,  USGS  Station  Number  08090800).  Wind 
speed  and  direction  are  routinely  measured  at  the  Granbury  Regional 
Airport  with  data  made  available  through  the  Local  Weather  Exchange 
Network  (station  KGDJ).  For  purposes  of  this  research,  wind  data  were 
converted  into  north-south  and  east- west  vectors. 

Project  components 

1.  Directly  measuring  water  exchange  rates  between  coves  and  natural 
surface  waters  of  the  open  lake,  and  exploring  relationships  between 
water  exchange,  wind  speed  and  direction,  and  inflow  to  the  lake. 
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Studies  were  conducted  in  three  coves,  where  the  water  exchange  rate  was 
measured  daily  over  a  period  of  5-7  days.  Each  cove  was  investigated 
thrice:  during  the  typical  periods  of  bloom  development  (February);  bloom 
decline  (April);  and  post-bloom  (June).  Daily  water  exchange  rates  were 
monitored  over  a  5-7-day  period  because  previous  research  employed  a 
similar  periodicity  (Hayden  et  al.  2012,  Roelke  et  al.  2012a). 

Study  coves  for  this  project  component  were  Ranger’s  Slough,  Bee  Creek, 
and  Thorp’s  Spring  (Figure  1).  These  coves  were  selected  based  on  their 
restricted  connection  to  the  mainstem  of  the  lake  and  their  larger  size. 
Their  restricted  connectivity  likely  increases  the  probability  of  longer 
hydraulic  residence  times,  thereby  increasing  the  probability  for  success  of 
future  bloom  mitigation.  Their  larger  size  likely  makes  for  more  effective 
refuge  habitat  able  to  shelter  larger  populations,  thereby  accelerating  the 
recovery  rate  of  populations  impacted  by  blooms  in  other  areas  of  the  lake. 


Figure  1.  Lake  Granbury  (USA)  where  5-  to  7-day-long  Rhodamine  WT 
dissipation  studies  were  conducted  in  three  coves  during  February, 
April,  and  June  2012. 


The  horizontal  and  vertical  dissipation  of  a  fluorescent  tracer,  Rhodamine 
WT  (Abbey  Color,  Inc.),  was  measured  in  order  to  estimate  the  water 
exchange  rate.  Rhodamine  WT  is  routinely  used  for  water  tracer  studies, 
and  is  benign  to  aquatic  organisms.  Because  of  its  fluorescing 
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characteristics,  it  was  detectable  by  fluorometer  at  very  low  concentrations 
(~0-5  ppb).  The  amount  of  Rhodamine  WT  needed  to  start  each  tracer 
experiment  was  determined  based  on  the  bathymetry  and  depth  of  the 
target  coves,  which  enabled  calculation  of  cove  volumes  (see  Dataflow 
methods  below).  To  start  each  experiment,  background  fluorescence  was 
first  determined  followed  by  tracer  addition.  The  tracer  was  added  by  slowly 
pumping  the  solution  into  the  boat  motor’s  outwash  as  the  boat  followed 
closely  spaced  transects  through  the  target  coves  (see  Figure  Ai).  This 
resulted  in  a  fairly  uniform  starting  condition. 

Horizontal  Rhodamine  WT  distribution  was  estimated  daily  using  high- 
resolution  spatial  mapping  technology.  Dataflow,  a  high-speed,  flow¬ 
through  measurement  apparatus  developed  for  mapping  physicochemical 
parameters  in  shallow  aquatic  systems  was  used  for  this  purpose  (Madden 
and  Day  1992).  The  Dataflow  configuration  allowed  concurrent  measure¬ 
ment  of  depth,  multiple  water  quality  parameters,  and  Rhodamine  WT. 
Measurements  were  taken  at  3-second  intervals  from  ~20  cm  below  the 
surface.  An  integrated  GPS  was  used  to  simultaneously  record  sample 
locations.  GPS  data  and  Dataflow  information  were  then  used  to  create 
gridded  data  files  (Surfer  v8.o).  These  data  files  were  then  used  to  create 
detailed  contour  maps  showing  the  surface  layer  distribution  of  Rhodamine 
WT  (Figures  A2-A10). 

Vertical  Rhodamine  WT  distribution  was  estimated  daily  by  vertical 
profiling  of  Rhodamine  WT  using  data  sondes  equipped  with  filter  sets 
designed  for  detection  of  this  tracer  (YSI,  Inc).  Multiparameter  data  sondes 
were  used  to  collect  basic  water  quality  information  (e.g.,  dissolved  oxygen, 
temperature,  pH,  specific  conductance),  which  allowed  for  determination  of 
chemoclines  and  thermoclines  to  support  numerical  modeling  of  coves. 
Fixed  station  locations  were  established  throughout  each  cove  for  this 
purpose  (see  Figure  Ai).  Another  gridded  data  file  was  then  generated 
expressing  the  ratio  of  total  Rhodamine  WT  in  the  water  column  to  the 
surface  Rhodamine  WT  concentration  at  each  of  the  fixed  stations. 

Total  mass  of  Rhodemine  WT  in  the  target  coves  was  determined  by 
combining  the  gridded  high-resolution  spatial  mapping  of  surface  waters 
with  the  gridded  vertical  profiling  data, using  the  gridding  functions  of  a 
commercial  software  package  (Surfer  v8.o).  An  exponential  decay  function 
was  then  fit  through  the  times  series  data  of  total  Rhodamine  WT 
(KaleidaGraph  4.0),  where  the  exponent  coefficient  reflected  the  rate  of 
water  exchange. 
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Relationships  between  cove-open  lake  water  exchange,  wind  speed  and 
direction,  and  inflow  were  investigated  using  multiple  linear  regression 
(SPSS  V20.0).  For  this  analysis,  the  24-hr  water  exchange  rate  for  each 
cove  was  calculated  using: 


E 


T>rl 

f2_t±2A_\ 

Rd 


t 


(1) 


where 

E  =  water  exchange  rate  (d1) 

Rdt  =  total  mass  of  Rhodamine  WT  in  the  target  cove  (kg) 

Rdt+24  =  total  mass  of  Rhodamine  WT  24  hrs  later  (kg) 
t  =  time  (d) 

Wind  and  inflow  data  used  in  this  analysis  were  the  24-hr  averages  prior 
to  the  time  of  sampling,  with  wind  data  split  into  its  north-south  and  east- 
west  vectors.  One  multivariate  linear  regression  was  performed  for  each 
target  cove,  i.e.  all  daily  estimates  of  water  exchange  rates  from  all  three 
experiments  were  gathered  into  one  analysis.  This  resulted  in  a  total  of  17, 
15,  and  12  observations  for  Rangers  Slough,  Bee  Creek,  and  Thorps  Spring, 
respectively. 

Note  that  this  approach  to  estimating  daily  water  exchange  assumed  that 
photolysis  and  absorption  of  Rhodamine  WT  over  24  hr  were  not  signifi¬ 
cant.  In  other  words,  reduction  in  total  Rhodamine  WT  in  the  coves  was  due 
to  hydraulic  dispersal  of  the  dye  into  the  open  lake  only.  Previous  studies 
quantifying  photolysis  and  absorption  of  Rhodamine  WT  concluded  that 
these  mechanisms  did  not  affect  estimation  of  hydraulic  processes  (Tai  and 
Rathbun  1988,  Upstill-Goddard  et  al.  2001,  Keefe  et  al.  2004,  Dierberg  and 
DeBusk  2005).  Regarding  photolysis,  it  was  concluded  that  accounting  for 
its  effects  would  be  necessary  if  experiments  lasted  several  months  to  years 
(Suijlen  and  Buyse  1994),  periods  much  longer  than  the  24-hr  period  used 
here. 

2.  Determining  seasonal  differences  between  plankton  community 
composition  and  water  quality  observed  between  coves  and  open  lake 
waters  in  Lake  Granbury  based  on  the  surmised  rate  of  natural  surface 
water  exchange. 
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The  study  used  historical  data  collected  in  Lake  Granbury  where  concurrent 
measurements  were  taken  at  to  stations  located  along  the  mainstem  of  the 
lake  (labeled  stations  l  through  10)  and  10  stations  located  in  connecting 
coves  (labeled  stations  A  through  J),  see  Figures  2  and  3.  Note  that  cove  G  is 
the  same  as  Ranger’s  Slough  and  cove  A  is  the  same  as  Bee  Creek.  Sampling 
occurred  monthly  during  this  previous  monitoring,  but  not  continuously, 
over  the  period  of  December  2007  through  May  2010.  Sampling  was  not 
performed  during  November  2008,  September  through  December  2009, 
and  March  2010. 

For  each  of  the  historical  sampling  trips,  parameters  measured  included 
characterizations  of  the  plankton  (total  phytoplankton  biomass  as 
chlorophyll  a,  P.  parvum  population  density  and  zooplankton  community 
composition),  nutrients  (nitrate  and  nitrite  summed,  NOx;  ammonium, 

NH4;  total  nitrogen,  TN;  soluble  reactive  phosphorus,  SRP;  and  total 
phosphorus,  TP)  and  water  quality  (temperature,  pH,  salinity,  and  Secchi 
depth).  This  historical  record  was  generated  from  multiple  previous  projects 
that  were  funded  by  Texas  Parks  and  Wildlife  (Brooks,  lead  PI)  and  the 
Department  of  Energy  (Roelke,  lead  PI). 


Figure  2.  Lake  Granbury  (USA)  where  monthly  monitoring  was  conducted 
from  December  2007  through  May  2010  at  fixed  location  stations  in  the 
open  lake  (filled  circles)  and  in  adjoining  coves  (open  squares).  Note, 
station  A  was  located  within  Bee  Creek  cove  and  station  G  was  located 
within  Ranger’s  Slough  cove. 
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Figure  3.  Depth  contours,  cove  volumes,  and  cove  mouth  cross-sectional  areas  for  the  coves 
of  Lake  Granbury  (USA)  from  which  historical  data  were  used.  These  locations  correspond  to 
sample  station  locations  identified  in  Figures  1  and  2  as  follows:  cove  A,  also  called  Bee 
Creek  (a);  cove  B  (b);  cove  C  (c);  cove  D  (d);  cove  E  (e);  and  cove  F  (f).  Locations  of  cove 
mouths  are  identified  with  blue  lines. 


To  complement  this  historical  data,  a  bathymetry  characterization  of  these 
coves  was  performed,  with  the  exception  of  cove  H  (due  to  an  orientation 
error).  This  was  achieved  using  the  Dataflow  apparatus  described 
previously.  Depth  data  collected  along  the  closely  spaced  boat  transects 
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through  each  of  these  coves  enabled  construction  of  depth  contour  maps. 
For  the  contour  maps,  cove  volumes  and  cross-sectional  areas  of  the  cove 
mouths  were  determined  (Surfer  v8.o). 


Figure  4.  Depth  contours,  cove  volumes,  and  cove  mouth  cross-sectional  areas  for  the  coves 
of  Lake  Granbury  (USA)  from  which  historical  data  were  used.  These  locations  correspond  to 
sample  station  locations  identified  in  Figures  1  and  2  as  follows:  cove  G,  also  called  Ranger’s 
Slough  (g);  cove  H  (h);  cove  I  (i);  and  cove  J  0)-  Cove  mouths  are  identified  with  blue  lines. 

The  Euclidean  distance  metric  was  used  to  estimate  the  dissimilarity 
between  coves  and  the  closest  station  location  from  the  mainstem  of  the 
lake.  This  allowed  researchers  to  quantitatively  assign  a  dissimilarity  value 
between  station  pairs  over  time,  where  greater  Euclidean  distances  signified 
greater  dissimilarity.  Prior  to  this  analysis,  data  were  transformed  by 
z-score. 

Data  from  cove  stations  A-E  were  paired  with  the  averaged  data  from  open 
lake  stations  2  and  3.  Data  from  cove  stations  F  and  G  were  paired  with  the 
averaged  data  from  open  lake  stations  4  and  5.  Data  from  cove  station  H 
were  paired  with  the  averaged  data  from  open  lake  stations  6  and  7;  data 
from  cove  station  I  were  paired  with  the  averaged  data  from  open  lake 
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stations  8  and  9;  and  data  from  cove  station  J  were  paired  with  data  from 
open  lake  station  10. 

Hierarchical  clustering  was  used  to  graphically  display  these  dissimilarities 
for  each  cove  and  the  nearest  open  lake  stations  (see  Appendix  C).  This 
resulted  in  10  cluster  analyses  (representing  each  cove-open  lake  station 
pair)  focused  on  the  biotic  parameters  and  10  cluster  analyses  focused  on 
the  abiotic  parameters.  Each  cluster  analysis  had  48  observations, 
representing  24  observations  in  the  cove  and  24  in  the  open  lake  over  the 
period  of  monitoring. 

In  the  resulting  hierarchy  of  clusters  represented  in  the  dendrograms,  the 
numbers  1-24  were  used  to  represent  the  time-series  sampling  of  the  cove 
station.  That  is,  T’  represented  the  sampling  in  the  cove  during  December 
2007,  ‘2’  represented  the  sampling  in  the  cove  during  January  2008,  and 
so  forth,  with  ‘24’  representing  the  last  cove  sampling  during  May  2010.  A 
similar  notation  was  followed  for  the  paired  open  lake  samplings  using  the 
numbers  25-48. 

The  dissimilarity  between  cove-open  lake  stations  was  then  compared  to 
wind  speed  and  direction,  and  inflow  to  the  lake.  This  relationship  was 
assessed  using  multiple  linear  regressions.  This  enabled  evaluation  of  the 
relative  importance  of  wind  speed  and  direction,  and  inflow.  This  analysis 
explored  all  combinations  of  winds  averaged  over  24,  48,  and  72  hr  prior 
to  the  time  of  sampling  and  inflows  averaged  over  7, 15,  and  30  days.  This 
document  reports  only  on  the  wind-inflow  combination  for  each  cove- 
open  lake  station  pair  that  explained  the  greatest  amount  of  variability. 
Finally,  scatter  plots  were  used  to  compare  the  strength  of  the  multivariate 
regression  models  to  the  cove  mouth  cross-sectional  area  to  volume  ratio. 

3.  Determining  the  sensitivity  of  P.  parvum  bloom  formation  under 

conditions  of  managed  deepwater  flushing  in  coves  of  Lake  Granbury  to 
the  rate  of  natural  surface  water  exchange  between  coves  and  the 
mainstem  of  Lake  Granbury. 

The  efficacy  of  deepwater  flushing  to  mitigate  P.  parvum  blooms  in  coves 
was  assessed  using  numerical  models  that  encompassed  the  range  of  water 
exchange  rates  documented  in  the  in-lake  tracer  experiment  studies 
(described  above).  These  modeling  studies  also  addressed  a  range  of 
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scenarios  for  other  biological  and  limnological  factors  that  affect  the 
population  dynamics  of  P.  parvum. 

The  starting  point  for  this  modeling  was  a  “chemostat”  (or  continuously- 
stirred  tank  reactor)  model  of  the  population  dynamics  of  P.  parvum  and 
varying  numbers  of  competing  algae.  This  model  permitted  forcing  of 
population  dynamics  by  known  relationships  of  P.  parvum  growth  rate  to 
temperature,  salinity,  and  nutrient  concentration  (Baker  et  al.  2007, 

2009).  This  model  was  previously  applied  to  population  dynamics  of  P. 
parvum  in  Lake  Granbury,  and  results  of  that  research  are  documented  in 
Grover  et  al.  (2010). Since  that  document  was  published,  modeling  has 
been  improved  by  adding  representations  of  toxin  production  by  P. 
parvum  and  zooplankton  grazing  on  algae  (Grover  et  al.  2011,  2012). 

For  this  study,  the  updated  model  was  embedded  in  the  “cove-main  lake” 
setting  developed  by  Grover  et  al.  (2011).  This  modeling  context  coupled 
two  chemostat  models,  one  representing  the  mainstem  of  a  large  lake  and 
the  other  a  cove,  and  allowed  for  setting  flows  in  both  segments  and 
exchanges  between  them,  to  represent  observed  rates.  Therefore  the  model 
used  is  a  coupled  biological -physical  model  focused  on  evaluating  bloom 
dynamics  within  coves. 

The  biological  component  of  the  model  included  several  interactions  and 
reactions  known  or  suspected  to  be  important  to  the  bloom  dynamics  of  P. 
parvum  (Figure  5).  Two  nutrients,  nitrate  and  phosphate,  are  consumed  by 
P.  parvum  and  other  algae,  and  support  their  growth.  The  other  algae, 
representing  cyanobacteria,  produce  a  dissolved  toxin  that  partially  inhibits 
the  growth  of  P.  parvum  (Roelke  et  al.  2010).  Both  types  of  algae  are 
consumed  by  a  small-bodied  zooplankton  population.  Consumption  of  the 
other  algae  supports  the  growth  of  the  zooplankton,  as  does  consumption  of 
P.  parvum,  but  the  dissolved  toxin  produced  by  P.  parvum  also  inhibits 
zooplankton  ingestion  and  growth.  The  relationship  between  zooplankton 
and  P.  parvum  therefore  ranges  from  positive  to  negative  depending  on  the 
amount  of  toxin  produced.  Toxin  production  by  P.  parvum  is  most  rapid 
when  its  population  growth  rate  is  reduced  by  suboptimal  temperature  or 
salinity  (Baker  et  al.  2007),  or  by  nutrient  limitation  (Johansson  and 
Graneli  1999,  Skingel  et  al.  2010).  Growth  rates  of  P.  parvum  and  other 
populations  depend  on  temperature,  salinity,  and  hydraulic  flow,  which  are 
provided  as  input  to  the  model  from  observational  data.  The  mathematical 
equations,  functions,  and  parameters  representing  these  relationships  are 
presented  in  detail  in  Appendix  B. 
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Figure  5.  Interactions  represented  by  the 
biological  reaction  components  of  the 
model.  Light  blue  arrows  indicate  positive, 
nutritional  interactions  and  red  arrows 
indicate  negative,  toxic  interactions. 

Equations  representing  the  reaction  rates  for  the  constituents  of  the 
biological  model  were  embedded  in  a  physical  model  representing  cove- 
main  lake  configuration  (Figure  6).  In  this  system,  a  cove  of  small  volume  is 
coupled  to  the  main  body  of  the  reservoir  by  an  exchange  flow  that  mixes 
constituents  between  these  compartments.  Throughflow  also  enters  and 
leaves  the  main  lake;  a  much  smaller  throughflow  enters  the  cove  and  exits 
into  the  main  lake.  The  biological  model  is  implemented  in  both  the  cove 
and  the  main  lake,  with  reaction  rates  that  depend  on  the  conditions  within 
each  respective  compartment.  An  additional  pumped  flow  is  imposed  in  the 
model  to  represent  hypothetical  bloom  treatments  that  enhance  the  small, 
natural  throughflow  of  the  cove  with  additional  flushing.  In  practice,  the 
source  water  for  this  flushing  would  be  deep  water  from  the  main  lake,  free 
of  P.  parvum,  which  tends  to  concentrate  in  the  upper  waters  of  the  lake. 

The  physical  component  of  the  model  is  a  set  of  ordinary  differential 
equations  representing  the  flows  illustrated  (Figure  6)  and  their  transport 
of  the  biological  constituents.  The  physical  parameters  quantifying  such 
transport  include  two  that  are  dimensionless  ratios:  ^is  the  ratio  of  cove 
volume  to  total  system  volume  (cove  plus  main  lake),  and  i//is  the  ratio  of 
flow  through  the  cove  to  total  flow  through  the  system.  Both  of  these  ratios 
are  small  (much  less  than  o.oi)  because  coves  are  small  in  relation  to  the 
main  lake.  For  simplicity,  these  ratios  are  set  as  equal.  While  this  is  not 
likely  to  be  strictly  true,  both  numbers  are  very  small,  and  as  long  as  this  is 
the  case,  differences  in  these  parameters  have  negligible  influence  on  the 
predictions  of  the  model.  Moreover,  the  detailed  hydrological  information 
necessary  to  estimate  i//is  not  available. 
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Figure  6.  Illustration  of  the  coupled  biological-physical  model  obtained 
by  embedding  the  biological  interactions  represented  in  Figure  5  into 
the  cove-main  lake  configuration.  Large,  dark  blue  arrows  indicate  flow 
through  the  main  lake,  and  smaller  dark  blue  arrows  indicate  exchange 
between  the  cove  and  the  main  lake  compartments.  The  dotted  dark 
blue  arrow  indicates  the  flow  that  is  imposed  by  pumping  deep  water 
from  the  main  lake  to  the  headwaters  of  the  cove. 


The  flow  through  the  total  system  (cove  plus  main  lake)  is  parameterized 
by  a  dilution  rate  D,  based  on  the  total  volume  and  total  throughflow  of  the 
system.  Daily  estimates  of  D  are  obtained  from  hydrological  data  on 
Brazos  River  flow  (Grover  et  al.  2010). 

An  additional  parameter  E  represents  the  exchange  rate  between  the  cove 
and  the  main  lake.  Estimation  of  this  parameter  was  the  goal  of  the  tracer 
experiments  done  in  this  study,  and  based  on  these  results,  a  plausible 
range  of  values  for  E  was  explored  in  modeling.  Although  E  varies  with  time, 
within  the  model  runs  reported  here,  its  value  was  held  constant  through 
time. 

Finally,  a  parameter  F  represents  the  additional  flushing  through  the  cove 
imposed  by  hypothetical  bloom  treatments. 

The  equations  of  the  physical  component  of  the  model  are: 
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Here,  the  variables  Xx  and  X2  are  the  concentrations  of  a  constituent  in  the 
main  lake  (compartment  l)  and  the  cove  (compartment  2),  respectively. 

The  first  term  in  each  equation  quantifies  the  inflow  of  a  constituent  into 
the  headwaters  of  the  main  lake  and  the  cove,  respectively.  For  the  main 
lake  and  the  cove,  it  is  assumed  phosphate  and  nitrate  enter  both  compart¬ 
ments  at  equal  concentrations  (i.e.  X"’  =  X2  )  and  values  for  these  inputs 

are  obtained  from  observations  during  P.  parvum  monitoring  studies 
(Grover  et  al.  2010,  2012).  Equal  input  concentrations  are  assumed  because 
fine-scaled  data  are  not  available  to  estimate  input  concentrations 
separately  for  coves. 

It  is  assumed  that  small  amounts  of  zooplankton,  P.  parvum,  and  other 
algae  are  input  into  the  main  lake  compartment  (i.e.  X"'  =  a  small  number) 

but  not  into  the  cove  (i.e.  X”1  =  0 ).  This  small  input  was  imposed  for 

technical  reasons,  to  avoid  extinction  of  any  species  during  long-term  model 
simulations  (Grover  et  al.  2010,  2012). 

Inputs  are  assumed  to  be  zero  for  both  the  toxin  produced  by  cyanobacteria 
and  the  toxin  produced  by  P.  parvum  (i.e.  X'"  =  X™  =  0  ).  It  was  also 

assumed  that  the  deep  water  used  for  flushing  has  zero  nutrient  concentra¬ 
tions.  This  is  unlikely  to  be  true,  but  because  the  model  does  not  completely 
budget  the  fate  of  the  nutrients  made  available  from  deep  water,  during  long 
model  runs  an  unrealistic  accumulation  occurs.  Thus,  for  essentially 
technical  reasons,  the  nutrient  content  of  deep  water  was  neglected  in 
modeling  reported  here. 

Prior  to  exploring  hypothetical  bloom  treatments  based  on  enhanced 
flushing  of  coves,  the  model  was  calibrated  against  observational  data  with 
the  flushing  parameter  E  set  to  zero.  The  calibration  procedure  was  similar 
to  that  used  for  the  one-compartment  lake  model  that  provided  the  basis  for 
the  biological  component  described  above  (Grover  et  al.  2012).  Calibration 
focused  on  a  one-year  time  period  beginning  1  August  2006.  Extensive 
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observations  of  P.  parvum  populations  and  other  lake  characteristics  are 
available  for  this  time  period.  Daily  values  of  temperature,  salinity,  light, 
nutrient  input,  and  dilution  rate  were  used  to  run  the  model  for  several 
(3-5)  annual  cycles,  producing  one-year  seasonal  cycles  of  the  predicted 
population  density  of  P.  parvum,  which  were  compared  to  observations. 

The  results  were  graphically  indistinguishable  for  a  range  of  parameters 
representing  the  volumes  and  exchange  rates  of  coves  as  measured  in  this 
study.  The  predicted  annual  cycle  of  P.  parvum  density  matched  that 
predicted  by  the  predecessor  model  (Grover  et  al.  2012),  and  thus  displays 
the  same  high  correspondence  to  observations  during  the  calibration  year 
(Figure  7).  Without  flushing  treatments,  the  dynamics  of  P.  parvum  in 
coves  are  predicted  to  be  virtually  identical  to  the  main  lake,  for  the  range 
of  exchange  rates  found  in  the  tracer  study.  That  is,  the  exchange  rates  are 
high  enough  to  predict  strong  coupling  of  cove  conditions  to  those  in  the 
main  lake. 


Day  of  the  Year 

Figure  7.  Dynamics  of  P.  parvum  in  the  main  lake  compartment  (solid  line)  predicted  by  the 
coupled  biological-physical  model  with  no  flushing  treatment  imposed,  after  calibration  to 

observations  (solid  circles). 

With  the  calibrated  model,  the  cove-specific  efficacy  of  mitigating  blooms 
with  flushing  treatments  was  evaluated  through  sensitivity  analyses  along 
a  gradient  of  water  exchange  rates  between  coves  and  the  mainstem  of  the 
lake,  where  the  gradient  explored  spans  the  level  of  water  exchange 
plausible,  based  on  the  tracer  studies  reported  here. 
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Model  runs  for  sensitivity  analysis  were  constructed  as  follows.  For  a  given 
cove,  the  model  was  run  with  no  flushing  treatment  for  four  years  using 
daily  values  of  temperature,  salinity,  light,  nutrient  input,  and  dilution  rate 
as  forcing  data.  This  was  a  sufficient  amount  of  time  to  achieve  the  annual 
cycle  illustrated  (Figure  7).  On  1  January  of  the  fifth  year  of  the  model  run, 
when  P.  parvum  density  exceeded  104  cells/ml,  flushing  treatments  began. 
That  is,  treatments  were  assumed  to  begin  after  establishment  of  a  large  P. 
parvum  bloom.  Treatments  continued  for  the  next  200  days,  which  covers 
the  early  part  of  the  year  when  P.  parvum  blooms  have  historically 
occurred  in  Lake  Granbury  and  other  Texas  reservoirs  (Roelke  et  al.  2011). 

The  design  of  the  sensitivity  analysis  was  based  on  plausible  scenarios  for 
flushing  treatments  during  the  200-day  treatment  period.  The  flows 
imposed  were  based  on  deploying  one,  two,  or  three  of  the  highest-capacity 
solar-powered  pumps  available  (SolarBee  model  SB10000HW  vi8, 
14.4-million  gallons/day  (mgd)  capacity).  Pumps  were  simulated  as 
operating  on  five  different  schedules:  operating  1,  2,  3,  5,  or  7  days  per  week. 
For  the  coves  that  were  studied  in  the  tracer  experiment,  deploying  this 
number  of  pumps  produces  flushing  that  equals  or  exceeds  the  flushing 
rates  found  to  be  effective  in  reducing  P.  parvum  blooms  in  enclosure 
experiments  (Hayden  et  al.  2012)  Intermittent  pumping  was  explored  in 
addition  to  continual  pumping  (operating  7  days  per  week)  because 
enclosure  experiments  imposed  flushing  treatments  intermittently. 
Therefore  a  total  of  15  bloom  treatment  scenarios  were  explored  for  each 
cove,  over  a  gradient  of  cove  exchange  rates  (E)  ranging  from  0.05  to  1.0  d1, 
going  both  above  and  below  the  exchange  rates  observed  in  tracer 
experiments. 

The  sensitivity  analysis  used  the  three  coves  in  which  tracer  experiments 
were  conducted.  Cove  volumes  measured  during  the  morphometric  survey 
prior  to  the  first  tracer  experiment  were  used  to  assign  the  volume  ratio 
parameter  (f>  (Table  1).  Although  cove  volumes,  and  hence,  values  of  (f>  can 
be  expected  to  vary  over  time  as  lake  levels  change,  this  variation  was 
neglected  because  the  resulting  changes  in  (f>  have  virtually  no  impact  on 
model  predictions.  The  flow  ratio  parameter  y/  was  set  equal  to  <f)  as 
explained  above.  The  flushing  parameter  F  was  calculated  from  the  pump 
capacity  (14.4  mgd  per  pump)  in  relation  to  cove  volume. 
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Table  1.  Volume  and  flow  characteristics  of  coves  simulated  for  the  sensitivity  analysis  of  the 

biological-physical  mode. 


Cove 

Volume  (m3) 

Volume  Ratio  (0) 

Flushing  (F,  d1)  based  on  number  of 
pumps 

1  Pump 

2  Pumps 

3  Pumps 

Bee  Creek 

140865 

0.000745 

0.3870 

0.7739 

1.1609 

Thorp’s  Spring 

243281 

0.00129 

0.2241 

0.4481 

0.6722 

Rangers  Slough 

357469 

0.00189 

0.1525 

0.3050 

0.4575 
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2  Technical  Reporting  and  Discussion 

Dye  tracer  studies 

The  slow  Rhodamine  WT  addition  method  adopted  in  this  study  resulted  in 
a  fairly  uniform  dye  distribution  throughout  the  target  coves.  A  representa¬ 
tive  example  is  shown  from  the  third  experiment  in  Ranger’s  Slough 
(Figure  8)  where  24  hr  after  dye  addition,  concentrations  of  Rhodamine  WT 
were  in  the  range  of  1300-1600  (pg)  throughout  most  of  the  cove.  The  daily 
horizontal  dissipation  of  dye  was  easily  followed  in  all  three  coves  during 
the  experiments  using  Dataflow  technology  (see  figures  in  Appendix  A). 
Similarly,  the  daily  vertical  dissipation  of  dye  was  easily  followed  from  the 
vertical  profiling  data.  Combined,  this  information  enabled  determination 
of  the  daily  dissipation  of  Rhodamine  WT  from  the  coves,  thereby  enabling 
estimation  of  water  exchange  rates  between  cove  and  open  lake  waters. 


32°27.00’ 
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Rhodamine  WT  (|_Lg) 


O 


32°26.l6’ 


-97°47-34’ 


-97°46.50’ 
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Figure  8.  Water  column  Rhodamine  WT  distribution  24  hr  after  addition  in  Ranger’s 
Slough,  Lake  Granbury  (USA)  during  the  third  experiment  (June  2012). 


Rangers  Slough  (cove  G,  from  below) 

Natural  surface  water  exchange  rates  during  the  three  experiments  were 
very  similar  in  Ranger’s  Slough,  with  the  exponent  coefficient  ranging 
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between  0.20  and  0.25  d1  (Figure  9).  Daily  variation  in  the  water  exchange 
rate  was  not  well  explained  by  24-hr  winds  and  lake  inflows.  The  regression 
model  showed  insignificant  p-values  for  these  independent  variables,  and 
insignificant  p-values  and  low  variance  are  explained  in  the  overall  model 
(Table  2). 


Figure  9.  Total  mass  of  Rhodamine  WT  in  Ranger’s 
Slough,  Lake  Granbury  (USA)  during  tracer  studies 
conducted  in  February  (a),  April  (b),  and  June  (c).  Flere, 
time-zero  marks  the  timing  of  Rhodamine  WT  addition 
immediately  after  measurement  of  background 
fluorescence. 
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Table  2.  Multiple  linear  regression  analysis  of  the  relationship  between  daily  water  exchange  rate, 
24-hr  wind  speed  and  direction,  and  24-hr  inflow  to  Lake  Granbury  for  the  three  Ranger’s  Slough 

experiments  combined. 


Cove 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Ranger’s  Slough 

NS_24 

0.805 

-0.144 

0.207 

0.698 

EW_24 

-0.098 

0.864 

lnflow_24 

-0.086 

0.791 

Bee  Creek  (cove  A,  from  below) 

In  Bee  Creek,  exponent  coefficients  were  again  similar  during  the  three 
experiments,  ranging  between  0.24  ando.30  d1  (Figure  3).  Again,  daily 
variation  in  the  water  exchange  rate  was  not  well  explained  by  24-hr  winds 
and  lake  inflows.  The  regression  model  showed  insignificant  p-values  for 
winds  and  inflow  variables,  and  insignificant  p-values  and  low  variance 
were  explained  in  the  overall  model  (Table  3). 

Thorp’s  Spring 

Natural  surface  water  exchange  rates  were  much  more  variable  during  the 
three  experiments  in  Thorp’s  Spring,  with  the  exponent  coefficient  ranging 
between  0.06  and  0.59  d1  (Figure  11).  The  low  rate  of  water  exchange 
during  the  third  experiment  might  have  been  an  artifact  of  a  passing  storm 
that  occurred  on  day  5  of  this  tracer  study.  The  storm  re-suspended  sedi¬ 
ments  that  likely  elevated  the  background  fluorescence  level,  resulting  in  an 
overestimation  of  Rhodamine  WT  late  in  the  experiment. 

Consistent  with  the  other  two  coves,  daily  variation  in  the  water  exchange 
rate  could  not  be  explained  by  24-hr  winds  and  lake  inflows.  The  regression 
model  showed  insignificant  p-values  for  these  independent  variables,  and 
lack  of  an  explanation  for  insignificant  p-values  and  low  variance  in  the 
overall  model  (Table  4). 

Analysis  of  historical  data 

Cove  A  (Bee  Creek,  from  above) 

Historical  data  were  used  to  analyze  the  dissimilarity  measure  between 
cove  station  A  and  the  average  of  open  lake  stations  2  and  3.  From  the 
biotic  parameters  (Figure  12)  it  can  be  seen  that  trends  in  biomass  over 
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Figure  10.  Total  mass  of  Rhodamine  WT  in  Bee  Creek, 
Lake  Granbury  (USA)  during  tracer  studies  conducted  in 
February  (a),  April  (b),  and  June  (c).  Here,  time-zero 
marks  the  timing  of  Rhodamine  WT  addition 
immediately  after  measurement  of  background 
fluorescence. 


Table  3.  Multiple  linear  regression  analysis  of  the  relationship  between  daily  water  exchange  rate, 
24-hr  wind  speed  and  direction,  and  24-hr  inflow  to  Lake  Granbury  for  the  three  Bee  Creek 

experiments  combined. 


Cove 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Bee  Creek 

NS_24 

0.597 

-0.080 

0.661 

0.191 

EW_24 

0.555 

0.283 

lnflow_24 

0.035 

0.917 
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Figure  11.  Total  mass  of  Rhodamine  WT  in  Thorp’s 
Spring,  Lake  Granbury  (USA)  during  tracer  studies 
conducted  in  February  (a),  April  (b),  and  June  (c).  Here 
time-zero  marks  the  timing  of  Rhodamine  WT  addition 
immediately  after  measurement  of  background 
fluorescence. 


Table  4.  Multiple  linear  regression  analysis  of  the  relationship  between  daily  water  exchange  rate, 
24-hr  wind  speed  and  direction,  and  24-hr  inflow  to  Lake  Granbury  for  the  three  Thorp’s  Spring 

experiments  combined. 


Cove 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Thorp’s  Spring 

NS_24 

0.671 

-0.145 

0.762 

0.270 

EW_24 

0.825 

0.251 

lnflow_24 

0.183 

0.613 
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Figure  12.  Biotic  parameters  from  the  historical  data  used  in  the  dissimilarity  analysis 
comparing  cove  A  (solid  lines)  with  the  averages  from  stations  2  and  3  (dashed  lines) 
including  P.  parvum  (a),  chlorophyll  a  (b),  heterotrophic  protozoa  (c),  copepod  adults  (d), 
copepod  nauplii  (e),  cladocera  (f),  and  rotifers  (g). 

time  were  similar  regarding  the  timing  in  peaks,  with  one-month  offsets 
sometimes  occurring.  Greater  variability  was  observed  in  the  magnitude  of 
peak  biomass  for  this  cove-open  lake  station  pair.  Abiotic  parameters 
(Figure  13)  showed  less  variability  for  this  cove-open  lake  station  pair, 
where  the  timing  and  magnitude  of  peaks  were  similar. 

As  stated  previously,  much  of  this  historical  data  was  generated  from 
projects  that  were  funded  by  Texas  Parks  and  Wildlife  (Brooks,  lead  PI)  and 
the  Department  of  Energy  (Roelke,  lead  PI)  and  has  been  reported 
previously.  These  data  have  been  included  in  this  section  only  to  familiarize 
the  reader  with  the  scope  of  data  analysis  being  performed  here. 
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Figure  13.  Abiotic  parameters  from  the  historical  data  used  in  the  dissimilarity  analysis 
comparing  cove  A  (solid  lines)  with  the  averages  from  stations  2  and  3  (dashed  lines) 
including  temperature  (a),  salinity  (b),  pH  (c),  Secchi  depth  (d),  nitrate  and  nitrite  summed, 
NOx  (e),  ammonium,  NH4  (f),  total  nitrogen,  TN  (g),  soluble  reactive  phosphorus  (h),  and  total 

phosphorus,  TP  (i). 
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Dendrograms  based  on  the  biotic  and  abiotic  parameters  showed  hierarchal 
clustering  where  the  dissimilarity  between  cove  and  open  lake  station  pairs 
could  be  evaluated  (Figure  14).  For  example,  the  dissimilarity  between 
biotic  parameters  of  the  cove  and  open  lake  during  December  2007  was 
~2  scaled  as  the  Euclidean  distance  (graphically  displayed  from  the  level  of 
clustering  between  observations  ‘1’  and  ‘25’).  The  dissimilarity  between  the 
cove  and  open  lake  changed  over  time,  and  did  not  show  consistent  seasonal 
patterns  (Figure  15). 

To  compare  cove  station  A  with  averaged  data  from  open-lake  stations 
2  and  3,  the  change  in  dissimilarity  (for  both  biotic  and  abiotic  parameters) 
over  the  monitoring  period,  a  sufficient  explanation  was  not  provided  by  the 
multiple  linear  regression  models.  None  of  the  wind  and  inflow  variables 
correlated  with  dissimilarity,  contributing  to  weak  overall  models  (Table  5). 


Figure  14.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  A  and  the 
averaged  data  from  stations  2  and  3.  Here,  red  numbers  reflect  observations  from  the  cove 
and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired 
based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open 
lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake 
sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 
open  lake  sampling  during  May  2010. 
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Figure  15.  Time-series  dissimilarity  between  biotic  (a) 
and  abiotic  (b)  parameters  as  estimated  by  the 
Euclidean  distance  for  the  comparison  between  cove 
station  A  and  open-lake  stations  2  and  3  averaged. 


Table  5.  The  best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and 
biotic  Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  A. 


Cove  A 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_48 

0.204 

0.081 

-0.041 

0.897 

EW_48 

-0.207 

0.513 

lnflow_15 

0.429 

0.045 

Biotic 

NS_24 

0.672 

-0.070 

-0.159 

0.572 

EW_24 

-0.253 

0.357 

lnflow_15 

0.170 

0.480 

Cove  B 

When  data  for  cove  station  B  were  compared  with  averaged  data  from 
open-lake  stations  2  and  3,  none  of  the  wind  and  inflow  variables 
correlated  well  with  dissimilarity  for  either  biotic  or  abiotic  parameters, 
and  the  multivariate  regression  models  were  weak  (Table  6). 

Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 
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Table  6.  Best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and  biotic 
Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  B. 


Cove  B 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_48 

0.191 

0.098 

-0.169 

0.615 

EW_48 

-0.430 

0.209 

lnflow_15 

0.389 

0.079 

Biotic 

NS_72 

0.571 

-0.047 

0.471 

0.200 

EW_72 

0.239 

0.494 

lnflow_15 

0.166 

0.493 

Cove  C 

For  the  biotic  parameters,  winds  along  the  east-west  vector  correlated  with 
dissimilarity  between  cove  station  C  and  averaged  data  from  open-lake 
stations  2  and  3.  This  resulted  in  a  stronger  multivariate  regression  model, 
although  the  p-value  still  exceeded  0.05  (Table  7). 


Table  7.  Best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and  biotic 
Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  C. 


Cove  C 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_72 

0.258 

0.059 

0.067 

0.827 

EW_72 

-0.436 

0.221 

lnflow_30 

-0.280 

0.277 

Biotic 

NS_24 

0.071 

0.193 

0.328 

0.232 

EW_24 

0.633 

0.020 

lnflow_30 

0.154 

0.474 

None  of  the  wind  and  inflow  variables  correlated  well  with  dissimilarity 
when  considering  the  abiotic  parameters,  and  the  multivariate  regression 
model  was  weak. 


Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 
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Cove  D 

Winds  along  the  north-south  vector  correlated  with  dissimilarity  between 
cove  station  D  and  averaged  data  from  open-lake  stations  2  and  3  when 
considering  biotic  parameters.  This  resulted  in  a  stronger  multivariate 
regression  model,  but  again  with  ap-value  still  exceeding  0.05  (Table  8). 


Table  8.  The  best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and 
biotic  Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  D. 


Cove  D 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_72 

<0.000 

0.609 

0.492 

0.017 

EW_72 

0.306 

0.159 

lnflow_30 

0.799 

<0.000 

Biotic 

NS_48 

0.150 

0.113 

0.702 

0.035 

EW_48 

0.409 

0.195 

lnflow_15 

0.220 

0.287 

Winds  along  the  north-south  vector  and  inflows  to  the  lake  correlated  with 
dissimilarity  when  considering  abiotic  parameters.  In  addition,  the  multiva¬ 
riate  regression  model  was  strong,  with  ap-value  less  than  0.05  and  ~6i% 
of  the  variance  explained. 

Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 

Cove  E 

For  the  biotic  parameters,  inflows  to  the  lake  correlated  with  dissimilarity 
between  cove  station  E  and  averaged  data  from  open-lake  stations  2  and  3. 
This  resulted  in  a  stronger  multivariate  regression  model,  although  with  a 
p-value  still  exceeding  0.05  (Table  9). 

None  of  the  wind  and  inflow  variables  correlated  well  with  dissimilarity 
when  considering  the  abiotic  parameters  (although  north-south  winds 
were  close),  and  the  multivariate  regression  model  was  insignificant. 

Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 
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Table  9.  The  best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and 
biotic  Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  E. 


Cove  E 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_72 

0.167 

0.102 

0.455 

0.076 

EW_72 

0.121 

0.624 

lnflow_15 

0.319 

0.133 

Biotic 

NS_72 

0.066 

0.199 

0.355 

0.140 

EW_72 

0.287 

0.274 

lnflow_30 

0.494 

0.040 

Cove  F 

North-south  and  east-west  wind  vectors,  as  well  as  inflows  to  the  lake, 
correlated  with  dissimilarity  between  cove  station  F  and  averaged  data 
from  open-lake  stations  4  and  5  when  considering  biotic  parameters.  A 
strong  multivariate  regression  model  resulted,  withp-value  less  than  0.05. 
However,  this  model  only  accounted  for  -36%  of  the  total  variability  in 
dissimilarity  (Table  10). 


Table  10.  The  best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and 
biotic  Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  F. 


Cove  F 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_72 

0.137 

0.128 

-0.295 

0.328 

EW_72 

-0.712 

0.031 

lnflow_15 

-0.195 

0.382 

Biotic 

NS_72 

0.009 

0.361 

2.798 

0.011 

EW_72 

2.113 

0.048 

lnflow_30 

2.600 

0.018 

For  the  abiotic  parameters,  the  east-west  wind  vector  correlated  with 
dissimilarity.  This  resulted  in  a  stronger  multivariate  regression  model, 
but  with  ap-value  still  exceeding  0.05. 


Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 
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Cove  G  (Ranger’s  Slough,  from  above) 

The  north-south  wind  vector  and  inflows  to  the  lake  correlated  with 
dissimilarity  between  cove  station  G  and  averaged  data  from  open-lake 
stations  4  and  5  when  considering  biotic  parameters.  A  strong  multivariate 
regression  model  resulted,  with  ap-value  less  than  0.05.  This  model 
accounted  for  -38%  of  the  total  variability  in  dissimilarity  (Table  11). 


Table  11.  The  best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and 
biotic  Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  G. 


Cove  G 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_48 

0.447 

-0.010 

0.239 

0.485 

EW_48 

0.433 

0.214 

lnflow_7 

0.246 

0.271 

Biotic 

NS_72 

0.009 

0.380 

0.673 

0.011 

EW_72 

0.430 

0.101 

lnflow_15 

0.452 

0.027 

None  of  the  wind  and  inflow  variables  correlated  well  with  dissimilarity, 
when  the  abiotic  parameters  and  the  multivariate  regression  model  were 
considered  to  be  insignificant. 

Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 

Cove  H 

The  east-west  wind  vector  and  inflows  to  the  lake  correlated  with 
dissimilarity  between  cove  station  H  and  averaged  data  from  open-lake 
stations  6  and  7  when  considering  biotic  parameters.  A  strong  multivariate 
regression  model  resulted,  with  ap-value  less  than  0.05.  This  model 
accounted  for  -58%  of  the  total  variability  in  dissimilarity  (Table  12). 

None  of  the  wind  and  inflow  variables  correlated  well  with  dissimilarity 
when  considering  the  abiotic  parameters  (although  inflows  to  the  lake 
were  close),  and  the  multivariate  regression  model  was  insignificant. 
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Table  12.  The  best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and 
biotic  Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  H. 


Cove  H 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_48 

0.198 

0.089 

0.228 

0.482 

EW_48 

0.507 

0.137 

lnflow_15 

0.419 

0.068 

Biotic 

NS_24 

<0.001 

0.576 

0.208 

0.287 

EW_24 

0.488 

0.018 

lnflow_15 

0.765 

0.000 

Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 

Cove  I 

When  considering  biotic  parameters,  inflows  to  the  lake  correlated  with 
dissimilarity  between  cove  station  I  and  averaged  data  from  open-lake 
stations  8  and  9.  A  strong  multivariate  regression  model  resulted,  with  a 
p-value  less  than  0.05.  This  model  accounted  for  -62%  of  the  total 
variability  in  dissimilarity  (Table  13). 


Table  13.  The  best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and 
biotic  Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  I. 


Cove  1 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_24 

0.840 

-0.104 

-0.099 

0.754 

EW_24 

-0.125 

0.697 

lnflow_15 

0.201 

0.383 

Biotic 

NS_24 

<0.000 

0.619 

0.066 

0.739 

EW_24 

0.297 

0.146 

lnflow_15 

0.827 

<0.000 

None  of  the  wind  and  inflow  variables  correlated  well  with  dissimilarity 
when  the  abiotic  parameters  and  the  multivariate  regression  model  were 
considered  to  be  insignificant. 
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Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 

Cove  J 

The  east-west  wind  vector  and  inflows  to  the  lake  correlated  with 
dissimilarity  between  cove  station  J  and  data  from  open-lake 
station  to  when  considering  biotic  parameters.  A  strong  multivariate 
regression  model  resulted,  with  ap-value  less  than  0.05.  This  model 
accounted  for  -41%  of  the  total  variability  in  dissimilarity  (Table  14). 


Table  14.  The  best  fit  multiple  linear  regression  analysis  of  the  relationship  between  the  abiotic  and 
biotic  Euclidean  dissimilarity,  wind  speed  and  direction,  and  inflow  to  Lake  Granbury  for  cove  J. 


Cove  J 

Predictor 

Variables 

Significance 
of  Model 

Adjusted 

R  square 

Standardized 

Coefficients 

Significance 
of  Predictor 
Variables 

Abiotic 

NS_72 

0.768 

-0.102 

-0.124 

0.770 

EW_72 

-0.272 

0.514 

lnflow_7 

0.151 

0.552 

Biotic 

NS_24 

0.008 

0.406 

0.635 

0.023 

EW_24 

0.672 

0.018 

lnflow_15 

0.557 

0.006 

None  of  the  wind  and  inflow  variables  correlated  well  with  dissimilarity 
when  the  abiotic  parameters  and  the  multivariate  regression  model  were 
considered  to  be  insignificant. 

Dendrograms  and  time-series  dissimilarity  graphs  corresponding  to  this 
cove-open  lake  pair  are  provided  in  Appendix  C. 

Role  of  cove  morphology 

When  comparing  the  strength  of  the  multivariate  regression  models  against 
the  morphology  of  coves,  correlations  were  poor  for  both  biotic  and  abiotic 
parameters  (Figure  16).  Here,  cove  morphology  is  only  depicted  by  the  cove 
mouth  cross-sectional  area  to  volume  ratio.  In  addition,  the  orientation  of 
coves  on  the  landscape  is  not  considered  here. 
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Figure  16.  Scatter  plots  with  linear  model  fits  comparing 
the  mouth  cross-sectional  area  to  volume  ratio  of  each 
cove  against  the  strength  of  the  multivariate  regression 
model  that  compared  the  dissimilarity  between  biotic  (a) 
and  abiotic  (b)  parameters  between  coves  and  open-lake 
stations  with  winds  and  inflows  to  the  lake. 


Numerical  modeling 

Rangers  Slough  (cove  G,  from  above) 

As  described  previously,  hypothetical  flushing  treatments  were  simulated 
for  a  range  of  cove  exchange  rates  encompassing  those  measured  in  the 
tracer  experiments,  for  each  of  three  coves.  A  coupled  biological-physical 
model  was  used,  treating  the  cove  and  the  main  lake  as  two  compartments 
coupled  by  hydraulic  exchange.  For  each  scenario  and  each  exchange  rate 
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assigned,  the  model  was  run  for  four  years  before  initiating  simulated 
treatments  on  January  l  of  the  fifth  year.  The  200  days  of  flushing  treat¬ 
ment  simulated  thereafter  followed  a  characteristic  pattern,  illustrated  here 
for  a  simulation  based  on  Rangers  Slough,  with  three  pumps  deployed,  and 
operated  for  either  5  or  7  days  each  week  (Figure  17).  For  this  illustration, 
the  lowest  measured  exchange  rate  (0.2  d1)  was  adopted. 


Rangers  Slough,  3  pumps,  5-7  days  on,  low  exchange 


Figure  17.  Illustrative  time  series  of  predicted  P.  parvum  population  density  in  Rangers 
Slough,  with  flushing  treatments  imposed  intermittently  (5  days  out  of  7),  or  continually 

(7  days  out  of  7). 


In  Figure  17,  three  time  series  of  predicted  P.  parvum  densities  are 
displayed.  The  upper  line  is  the  predicted  population  density  in  the  main 
lake.  The  other  two  lines  display  predicted  population  densities  in  the 
cove.  Time  is  measured  in  days  from  the  start  of  the  fourth  year  of  the  run, 
and  all  three  lines  coincide  until  day  365,  the  start  of  the  fifth  year,  when 
flushing  treatments  are  imposed.  After  the  lines  diverge,  the  lower  line  is 
the  predicted  population  density  with  continual  flushing  of  the  cove 
(7  days  of  pumping  per  week).  The  intermediate  line,  bouncing  from  high 
to  low,  is  the  predicted  population  density  with  intermittent  flushing  of  the 
cove  (5  days  of  pumping  per  week).  The  model  predicts  that  when  pumps 
are  on,  the  population  density  of  P.  parvum  is  rapidly  reduced  to  a  lower 
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value.  However,  with  intermittent  pumping,  population  density  recovers 
quickly  during  the  two  days  that  pumps  are  off  and  is  close  to  that  in  the 
main  lake. 

These  results  indicate  that  continual  flushing  is  predicted  to  reduce  the 
density  of  P.  parvum.  However,  intermittent  flushing  is  rapidly  overcome 
by  exchange  between  the  cove  and  main  lake,  such  that  P.  parvum  is 
reduced  when  pumping  is  on,  but  rapidly  rises  when  it  is  off.  These  large 
fluctuations  are  driven  by  the  force  of  exchange  with  the  much  larger  main 
lake,  and  arise  even  when  a  relatively  modest  exchange  rate  is  assumed.  To 
make  long-term  comparisons  among  treatment  scenarios  using  different 
numbers  of  pumps  and  schedules  of  intermittent  (or  continual)  operation, 
P.  parvum  density  over  the  200  days  of  simulated  flushing  treatment  was 
averaged,  and  the  percent  reduction  below  the  main  lake  level  was 
calculated  for  each  model  run  and  displayed  in  relation  to  the  exchange 
rate  (Figures  18  and  19). 

The  predicted  reduction  of  P.  parvum  is  sensitive  to  the  treatment  scenario 
adopted,  and  to  the  value  of  the  exchange  rate  (Figures  18  and  19).  For  any 
scenario,  the  reduction  is  predicted  to  decrease  as  the  exchange  rate 
increases.  This  is  expected,  because  with  a  larger  exchange  rate,  pumping 
must  overcome  more  rapid  transport  of  P.  parvum  from  the  main  lake.  For 
the  observed  range  of  exchange  rates  in  Rangers  Slough,  0.20  to  0.25  d-i, 
reduction  in  P.  parvum  relative  to  the  main  lake  ranges  from  about  5%  for 
the  least  effective  treatments,  up  to  about  70%  for  the  most  effective. 
Increasing  the  number  of  pumps  or  the  number  of  days  that  they  run  always 
produces  a  greater  reduction  of  P.  parvum.  For  scenarios  using  1,  2,  or  3 
pumps,  continual  pumping  is  always  more  effective,  and  reduces  P.  parvum 
density  by  40-70%  at  the  range  of  exchange  rates  observed  in  Rangers 
Slough. 

Bee  Creek  (cove  A,  from  above) 

The  results  of  modeling  various  bloom  treatment  scenarios  for  Bee  Creek 
are  similar  to  those  described  previously  for  Rangers  Slough.  The  predicted 
reduction  of  P.  parvum  is  again  sensitive  to  the  treatment  scenario  adopted, 
and  to  the  value  of  the  exchange  rate  (Figures  20  and  21).  For  any  scenario, 
the  reduction  is  predicted  to  decrease  as  the  exchange  rate  increases.  For 
the  observed  range  of  exchange  rates  in  Bee  Creek,  0.24  to  0.30  d1, 
reduction  in  P.  parvum  relative  to  the  main  lake  ranges  from  about  8%  for 
the  least  effective  treatments,  up  to  about  80%  for  the  most  effective. 
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Rangers  Slough  -  one  pump  (14.4  mgd) 


Rangers  Slough  -  two  pumps  (28.8  mgd) 


Figure  18.  Reduction  of  P.  parvum  in  Rangers  Slough 
for  flushing  treatment  scenarios  with  intermittent 
pumping  for  short  periods  (1-3  days  out  of  7). 
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Rangers  Slough  -  one  pump  (14.4  mgd) 


Rangers  Slough  -  two  pumps  (28.8  mgd) 


Exchange  between  main  lake  and  cove  (per  day) 

Figure  19.  Reduction  of  P.  parvum  in  Rangers  Slough 
for  flushing  treatment  scenarios  with  intermittent 
pumping  for  long  periods  (3-5  days  on  out  of  7),  or 
continual  pumping  (7  days  on  out  of  7). 
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Bee  Creek  -  one  pump  (14.4  mgd) 


Bee  Creek  -  two  pumps  (28.8  mgd) 


Figure  20.  Reduction  of  P.  parvum  in  Bee  Creek  for 
flushing  treatment  scenarios  with  intermittent  pumping 
for  short  periods  (1-3  days  out  of  7). 
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Bee  Creek  -  one  pump  (14.4  mgd) 


Exchange  between  main  lake  and  cove  (per  day) 

Figure  21.  Reduction  of  P.  parvum  in  Bee  Creek  for 
flushing  treatment  scenarios  with  intermittent  pumping 
for  long  periods  (3-5  days  out  of  7),  or  continual  pumping 
(7  days  out  of  7). 

Increasing  the  number  of  pumps  or  the  number  of  days  that  they  run  always 
produces  a  greater  reduction  of  P.  parvum.  For  scenarios  using  l,  2,  or  3 
pumps,  continual  pumping  is  always  more  effective,  and  reduces  P.  parvum 
density  by  60-80%  at  the  range  of  exchange  rates  observed  in  Bee  Creek. 
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Thorp’s  Spring 

The  results  of  modeling  various  bloom  treatment  scenarios  for  Thorp’s 
Spring  are  similar  to  those  described  previously  for  Rangers  Slough  and  Bee 
Creek.  The  predicted  reduction  of  P.  parvum  is  again  sensitive  to  the  treat¬ 
ment  scenario  adopted,  and  to  the  value  of  the  exchange  rate  (Figures  22 
and  23).  For  any  scenario,  the  reduction  is  predicted  to  decrease  as  the 
exchange  rate  increases.  For  the  observed  range  of  exchange  rates  in 
Thorp’s  Spring,  0.38  to  0.59  d1  (excluding  a  low  value  of  0.06  d1  whose 
measurement  is  suspect),  reduction  in  P.  parvum  relative  to  the  main  lake 
ranges  from  about  5%  for  the  least  effective  treatments,  up  to  about  65%  for 
the  most  effective.  Increasing  the  number  of  pumps  or  the  number  of  days 
that  they  run  always  produces  a  greater  reduction  of  P.  parvum.  For 
scenarios  using  1,  2,  or  3  pumps,  continual  pumping  is  always  more 
effective,  and  reduces  P.  parvum  density  by  30-65%  at  the  range  of 
exchange  rates  observed  in  Thorp’s  Spring. 
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Thorp  Spring  -  one  pump  (14.4  mgd) 


Cove  -  Main  Lake  Exchange  Rate  (per  day) 

Figure  22.  Reduction  of  P.  parvum  in  Thorp’s  Spring  for 
flushing  treatment  scenarios  with  intermittent  pumping 
for  short  periods  (1-3  days  out  of  7). 
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Thorp's  Spring  -  one  pump  (14.4  mgd) 


Thorp's  Spring  -  two  pumps  (28.8  mgd) 


|  100 


0.2  0.4  0.6  0.8  1.0 


Exchange  between  main  lake  and  cove  (per  day) 

Figure  23.  Reduction  of  P.  parvum  in  Thorp’s  Spring  for 
flushing  treatment  scenarios  with  intermittent  pumping 
for  long  periods  (3-5  days  out  of  7),  or  continual  pumping 
(7  days  out  of  7). 
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3  Discussion 

Efficacy  of  deepwater  flushing 

Natural  water  exchange  rates  between  coves  and  surface  waters  of  the  open 
lake  were  high.  For  Ranger’s  Slough  and  Bee  Creek,  two  coves  with  a  lower 
cove  mouth  cross-sectional  area  to  volume  ratio,  water  exchange  rates  were 
-0.25  d1.  For  Thorp’s  Spring,  exchange  rates  were  more  variable.  But  the 
one  low  exchange  rate  observed  for  this  cove  coincided  with  a  late- 
experiment  wind  event  that  suspended  sediments.  These  sediments  likely 
interfered  with  on-board  fluorometer  measurements.  In  turn,  this  likely 
resulted  in  an  overestimation  of  Rhodamine  WT  late  in  the  dye  tracer 
experiment,  thereby  masking  the  true  dissipation  rate.  If  that  observation  is 
discounted,  then  water  exchange  rates  for  Thorp’s  Spring  were  higher  than 
the  other  two  coves.  This  finding  would  be  consistent  with  having  a  larger 
cove  mouth  cross-sectional  area  to  volume  ratio  for  Thorp’s  Spring  (0.0014) 
compared  to  Ranger’s  Slough  (0.0007)  and  Bee  Creek  (0.0011). 

Given  water  exchange  rates  of  -0.25  d1  and  higher,  it  is  not  surprising 
that  plankton  dynamics  and  water  quality  changes  in  the  coves  of  Lake 
Granbury  were  similar  to  those  of  the  open  lake  in  the  historical  record 
(Roelke,  Brooks,  and  Grover,  unpublished  data).  Homogenization  effects 
in  regards  to  community  composition  occur  in  systems  where  immigration 
rates  are  high  (Leibold  and  Norberg  2004,  Leibold  et  al.  2004).  They  occur 
because  spatial  heterogeneity  in  the  system  is  lost.  In  other  words,  the 
system  behaves  as  a  single  closed  entity  driven  by  processes  acting  at  the 
system  scale,  which  overwhelm  local-scale  processes.  This  generalized 
system  response  to  high  immigration  has  been  observed  in  many  other 
ecosystem  types  that  include  sea  floor  communities,  alpine  vegetation,  and 
island  fauna  (Thrush  et  al.  2006,  Britton  et  al.  2009,  Massed  2009).  In 
Lake  Granbury,  it  is  likely  that  spatial  heterogeneity  is  reduced  due  to  high 
water  exchange  rates  between  coves  and  the  surface  waters  of  the  open 
lake.  The  resulting  high  immigration  of  plankton  likely  leads  to  similar 
seasonal  patterns  observed  between  coves  and  the  open  lake. 

In  addition  to  homogenization  effects  arising  at  these  high  water  exchange 
rates,  it  is  likely  that  the  presence  of  P.  parvum  influences  plankton 
community  composition  and  the  seasonal  plankton  trajectory.  This  alga 
directly  affects  many  competing  phytoplankton  through  allelopathy,  many 
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zooplankton  predators  through  toxicity,  and  preys  on  heterotrophic  bacteria 
and  some  small  protozoa  (Ulitzur  and  Shilo  1964,  Igarashi  et  al.  1995, 
Nygaard  and  Tobiesen  1993,  Tillmann  2003,  Fistarol  et  al.  2003,  Rosetta 
and  McManus  2003,  Barreiro  et  al.  2005,  Baker  et  al.  2007,  2009).  Some  of 
these  interactions  were  observed  at  sub-lethal  levels,  and  some  effects  were 
taxon-specific  (Roelke  et  al.  2007,  2012b;  Brooks  et  al.  2010,  Remmel  et  al. 
2011).  So  it  is  likely  that  P.  parvum  influences  plankton  compositional 
states,  reducing  the  probability  of  differing  community  compositions,  even 
during  non-bloom  periods.  In  this  sense,  P.  parvum  would  be  serving  in  the 
role  of  an  ecosystem-engineering  organism  ( sensu  Jones  et  al.  1994). 

Interestingly,  water  exchange  rates  between  coves  and  surface  waters  of 
the  open  lake  were  not  well  predicted  by  the  24-hr  averaged  wind  speed 
and  direction,  or  inflow  to  the  lake.  It  may  be  that  dissipation  of  turbu¬ 
lence  important  to  cove-open  lake  water  exchanges  occurs  over  longer 
periods.  Alternatively,  it  may  be  that  biotic  and  abiotic  attributes  of 
plankton  environments  are  sensitive  to  processes  occurring  over  periods 
less  than  24  hr,  such  as  grazing  and  excretion  events  by  zooplankton  and 
planktivorus  fish.  An  experimental  design  involving  several  more  dye- 
tracer  experiments  timed  throughout  the  year  during  periods  of  higher 
and  lower  winds  and  inflows  would  be  needed  to  explore  this  directly. 
Currently,  these  data  do  not  exist  for  Lake  Granbury.  Indirect  evidence, 
however,  showed  that  winds  and  inflows  over  longer  periods  influenced 
water  exchanges  between  some  coves  and  the  open  lake  in  the  analysis  of 
historical  data  (discussed  below). 

As  mentioned  previously,  plankton  dynamics  and  water  quality  changes  in 
the  coves  of  Lake  Granbury  were  similar  to  those  of  the  open  lake.  But  they 
were  not  identical.  The  timing  and  magnitude  of  biomass  and  chemical 
concentration  peaks  varied  with  temporal  offsets  of  ~i-2  months  and  peak 
magnitudes  differed  several-fold  at  times.  Dissimilarity  analysis  focuses  on 
these  types  of  differences. 

Temporal  changes  in  dissimilarity  measures  were  dynamic,  with  some  coves 
showing  greater  variability  than  others.  For  example,  cove  A  (Bee  Creek) 
showed  a  gradual  increase  in  dissimilarity  with  the  open  lake  in  2008  while 
cove  G  (Ranger’s  Slough)  showed  erratic  changes  in  dissimilarity  during  this 
same  period.  Intuitively,  and  in  reference  to  biotic  parameters,  one  can 
envision  dissimilarity  increasing  over  time  after  local  communities  become 
isolated,  a  notion  stemming  from  Island  Biogeography  Theory  (MacArthur 
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and  Wilson  1967).  In  plankton  systems,  where  reproduction  and 
competitive  exclusion  occur  over  short  periods,  it  is  anticipated  that 
dissimilarity  between  communities  might  arise  quickly,  even  at  aggregated 
taxonomic  levels  (Cottenie  et  al.  2003,  Xu  et  al.  2012).  With  the  onset  of 
immigration  that  would  arise  when  waters  begin  to  mix  again,  dissimilarity 
can  decrease.  This  phenomenon  was  suggested  to  happen  very  quickly  in 
plankton  systems,  and  at  very  low  levels  of  immigration  (Roelke  and 
Eldridge  2008,  2010).  Dissimilarity  would  all  but  disappear  with  water 
exchanges  between  the  cove  and  surface  waters  of  the  open  lake,  leading  to 
a  homogenized  system  state. 

It  is  expected  that  water  exchanges  will  be  driven  primarily  by  winds  and 
inflows  to  a  lake.  As  such,  it  is  anticipated  that  dissimilarity  between  coves 
and  the  open  lake  should  also  closely  follow  winds  and  inflows.  For  the 
shorter  periods  investigated  here  (24,  48,  and  72  hr),  this  was  only  observed 
for  some  coves,  however.  For  example,  coves  F,  G,  H,  I,  and  J,  all  coves 
lower  in  the  lake,  showed  a  significant  relationship  in  the  multivariate 
regression  models  of  biotic  parameters.  These  models  were  strongly  driven 
by  inflows  to  the  lake.  Coves  closer  to  the  head  of  the  lake  did  not  show  a 
significant  relationship  in  the  multivariate  regression  models. 

It  may  be  that  inflows  to  the  lake  follow  a  laminar  pattern  of  flow  initially, 
i.e.,  when  in  “upstream”  areas  of  the  lake,  but  eventually  yield  to  a  more 
turbulent  flow  regime  in  “downstream”  areas  of  the  lake.  Under  such  a 
scenario,  lateral  water  exchanges  with  adjacent  coves  would  be  enhanced 
lower  in  the  lake,  and  those  coves  would  be  more  vulnerable  to  homogeniza¬ 
tion  effects.  The  multivariate  models  that  use  the  biotic  data  fit  this  explana¬ 
tion.  The  models  using  abiotic  data,  however,  did  not  show  this  same  trend. 
This  is  not  surprising,  however,  in  that  abiotic  parameters  would  not  show 
sensitivity  to  mixing  on  par  with  the  expected  sensitivity  of  biotic 
parameters  to  immigration. 

Winds  also  showed  significant  relationships  in  the  multivariate  regression 
models  of  biotic  and  abiotic  parameters.  But  no  consistent  spatial  trend 
emerged.  For  example,  winds  blowing  in  the  direction  of  a  cove  mouth 
sometimes  were  significant  in  the  multivariate  regression  models,  as  in  cove 
G  (Ranger’s  Slough).  Surprisingly,  this  was  the  exception,  not  the  norm.  In 
addition,  winds  were  commonly  significant  in  the  multivariate  regression 
models  when  they  were  blowing  perpendicular  to  the  orientation  of  the 
cove’s  mouth.  It  is  possible  Ekman  transport  might  have  developed  in  such 


ER  DC/EL  CR-13-1 


47 


cases,  driving  a  water  exchange  between  cove  and  lake  surface  waters 
perpendicular  to  the  direction  of  winds.  Exploration  of  this  circulation 
mechanism  is  beyond  the  scope  of  this  study,  however. 

Dissimilarity  in  the  coves  situated  higher  in  the  lake  (coves  A-D)  was  not 
well  explained  in  the  multivariate  regression  models,  and  the  coves  highest 
in  the  lake  (coves  A  and  B)  never  correlated  with  winds  or  inflows.  Further¬ 
more,  the  strength  of  the  multivariate  regression  models  was  poorly  related 
to  the  measure  of  cove  morphometry,  i.e.,  cove  mouth  cross-sectional  area 
to  cove  volume  ratio. 

Other  factors  likely  influenced  water  exchange  between  coves  and  the  open 
lake.  These  include  morphometry  of  the  open  lake  adjacent  to  coves  and 
the  lake’s  variable  fetch.  In  regards  to  open  lake  morphometry,  sills  have 
developed  over  the  life  span  of  this  reservoir  and  may  impede  water 
movement  towards  some  coves.  Sills  may  be  more  likely  to  form  laterally 
to  the  mainstem  of  the  lake  closer  to  the  point  of  river  inflow.  In  that  area, 
current  deceleration  begins  as  the  system  widens  and  deepens,  allowing 
sediment  deposition.  In  regards  to  fetch,  Lake  Granbury  is  sinuous.  This 
results  in  coves  with  the  same  landscape  orientation  having  varied 
relationships  with  winds.  For  example,  cove  I  would  be  more  susceptible 
to  wind-driven  mixing  with  southerly  winds  relative  to  cove  C  because  of 
the  open  expanse  of  lake  stretching  southward  from  cove  I  and  the  narrow 
section  of  open  lake  in  front  of  cove  C.  It  does  not  seem  likely  that 
differences  in  water  density  would  have  an  effect,  as  water  temperatures 
were  very  similar  between  the  coves  and  open  lake,  even  when  considering 
deeper  waters,  at  this  time  of  year. 

While  the  results  of  this  study  shed  only  some  light  on  the  mechanisms 
driving  water  exchange  between  coves  and  surface  waters  of  the  open  lake, 
they  conclusively  show  that  measured  rates  of  water  exchange  are  high  and 
that  the  coves  of  Lake  Granbury  vary  in  their  morphometry.  Consequently, 
it  follows  that  mitigating  P.  parvum  blooms  in  coves  through  deepwater 
flushing  will  require  different  levels  of  pumping. 

The  model  enables  determination  of  pumping  volumes  and  the  duration 
pumps  need  to  be  turned  on  to  achieve  varied  levels  of  bloom  reduction.  For 
example,  at  the  0.25  d1  water  exchange  rate  and  assuming  a  target  50% 
reduction  in  P.  parvum  bloom,  the  model  estimates  that  ~29  mgd  of 
deepwater  flushing  (two  pumps)  for  5  days  per  week  would  be  required  in 
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Ranger’s  Slough.  Pumping  at  the  lower  rate  of  -15  mgd  (one  pump)  would 
not  attain  the  targeted  50%  bloom  reduction.  For  Bee  Creek,  a  50%  bloom 
reduction  could  be  achieved  by  pumping  -29  mgd  for  4  days  per  week,  or  by 
pumping  -15  mgd  for  5  days  per  week.  For  Thorp’s  Spring,  a  50%  bloom 
reduction  could  be  achieved  by  pumping  -29  mgd  for  5  days  per  week  or  by 
continuously  pumping  at  -15  mgd. 

The  measured  exchange  rates  between  coves  and  the  open  lake  are  large 
enough  so  that  main  lake  and  cove  conditions  are  coupled  in  the  absence 
of  pumping,  with  homogenization  effects  masking  the  influence  of 
localized  processes.  Modeling  suggests  that  the  coupling  is  not  so  strong 
that  it  cannot  be  overcome  by  deepwater  pumping  to  reduce  P.  parvum. 
The  model  also  suggests  that  pumping  scenarios  with  higher  and  more 
continual  flows  will  be  the  most  effective. 

Recommendations 

Results  of  this  study  suggest  that  P.  parvum  blooms  can  effectively  be 
mitigated  in  cove  habitats  of  Lake  Granbury  through  managed  deep-water 
hydraulic  flushing.  The  acting  mechanisms  that  enable  this  approach  to  be 
effective  are  the  hydraulic  displacement  of  P.  parvum  populations  from 
cove  waters  and  the  ecological  disruption  of  bloom  development  that  occurs 
with  nutrient  loading  associated  with  flushing  events.  The  next  steps  in 
research  aimed  at  further  developing  this  mitigation  approach  should  be  a 
coupled  study  involving  model  refinement  and  a  valuation  study. 

Model  refinement 

Model  predictions  probably  underestimate  the  flushing  magnitude  and 
duration  of  flow  needed  to  mitigate  P.  parvum  blooms.  They  assume  the 
coves  to  be  well-mixed  entities.  They  are  not.  A  model  capturing  the 
longitudinal  aspects  of  coves,  perhaps  employing  a  one-dimensional  frame¬ 
work,  would  better  represent  cove  hydrology.  Within  a  one-dimensional 
framework  for  the  coves,  the  required  amount  of  deep-water  flushing  from 
the  open  lake  would  be  less  than  the  estimates  reported  here.  When  a  cove 
is  flushed,  its  upper  waters  are  likely  cleared  of  P.  parvum  to  a  high  degree, 
even  if  the  mouth  of  the  cove  retains  higher  populations. 

In  addition,  the  open  lake  should  be  modeled  with  greater  detail.  While 
the  open  waters  are  not  stratified  in  regards  to  temperature  at  this  time  of 
year,  chemoclines  are  still  likely  to  exist.  Exchanging  water  between  the 
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open  lake  and  coves  at  the  volumes  suggested  here  might  have  an  effect  on 
water  chemistry  of  the  open  lake  and  possibly  influence  water  quality 
“downstream.”  For  this  purpose,  off-the-shelf  models  like  CE-QUAL-W2 
might  be  used. 

Valuation  study 

The  costs  of  infrastructure,  pump  operation,  and  maintenance  need  to  be 
determined.  The  cost  of  pipe  may  be  high,  as  it  will  have  to  contain  waters 
being  transported  from  deeper  waters  of  the  open  lake  to  the  head  regions 
of  coves.  Lower  pumping  levels  will  lead  to  smaller  diameter  pipe  require¬ 
ments,  and  smaller  pumps  with  lower  energy  consumption,  thereby 
lowering  total  costs.  This  underscores  the  importance  of  model  refinement. 
A  valuation  study  would  then  compare  these  costs  to  what  will  be  gained  by 
implementing  this  mitigation  strategy.  Here,  gains  include  an  accelerated 
recovery  rate  of  fisheries  lake-wide,  renewal  of  recreational  use  of  the  lake, 
and  boosted  commerce  associated  with  increased  visitors  to  the  region.  On 
these  latter  points,  recovery  rate  of  fisheries  lake-wide  and  renewal  of 
recreational  use  of  the  lake  following  successful  P.  parvum  mitigation  in 
coves  has  yet  to  be  studied. 

Need  for  an  in-lake  demonstration 

The  structure  and  parameterization  of  the  numerical  model  is  based  on  a 
decade  of  research  focused  on  P.  parvum  in  Texas  waters.  Even  so,  there  are 
many  uncertainties  regarding  ecological  and  foodweb  interactions  prior  to 
and  during  blooms.  For  example,  there  appear  to  be  grazers  that  are 
resistant  to  P.  parvum  toxins,  and  there  appear  to  be  bacterial  and/or  viral 
pathogens  that  affect  P.  parvum.  Neither  of  these  were  depicted  in  the 
modeling  due  to  lack  of  information,  and  either  may  be  important  as  a 
bloom-influencing  factor.  Continued  research  is  needed  to  focus  on  issues 
such  as  these,  as  it  will  improve  the  model  predictions.  But  all  models  are 
ultimately  limited  in  that  they  can  never  capture  all  the  complexities  of  a 
natural  ecosystem.  As  such,  models  can  only  predict  what  might  happen. 
Consequently,  an  in-lake  demonstration  is  needed  for  situations  where 
ecosystem  complexities,  even  if  unknown,  are  influencing  the  system.  Such 
a  demonstration,  guided  by  results  from  a  more  refined  model  and  valua¬ 
tion  study,  is  needed  to  advance  the  notion  of  P.  parvum  mitigation  in 
Texas  lakes. 
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Appendix  A:  High-resolution  Spatial  Mapping 
of  Rhodamine  WT 

This  appendix  shows  the  high-resolution  spatial  mapping  of  Rhodamine 
WT  in  the  three  coves  of  focus  for  the  tracer  studies:  Ranger’s  Slough,  Bee 
Creek,  and  Thorp’s  Spring.  Also  shown  are  the  transects  the  boat  followed 
while  generating  the  data  for  the  high-resolution  maps,  and  the  fixed 
position  stations  where  the  vertical  profiling  was  performed. 
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Figure  Al.  Rhodamine  WT  dissipation  experiments  were  conducted  in  Ranger’s  Slough  (a), 
Bee  Creek  (b),  and  Thorp’s  Spring  (c),  Lake  Granbury  (USA).  Experiments  consisted  of  surface 
water  horizontal  mapping  based  on  a  boat  following  closely  spaced  transects  (dotted  lined) 
and  vertical  profiling  at  fixed  location  stations  (filled  circles). 
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Figure  A2.  Rhodamine  WT  distribution  in  Ranger’s  Slough  during  the  first  experiment 
(February)  on  days  one  (a),  two  (b),  three  (c),  four  (d),  five  (e),  and  six  (f). 
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Figure  A3.  Rhodamine  WT  distribution  in  Ranger’s  Slough  during  the  second  experiment 
(April)  on  days  one  (a),  two  (b),  three  (c),  and  four  (d). 
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Figure  A4.  Rhodamine  WT  distribution  in  Ranger’s  Slough  during  the  second  experiment 
(April)  on  days  five  (e),  six  (f)  and  seven  (g). 
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Figure  A5.  Rhodamine  WT  distribution  in  Ranger’s  Slough  during  the  third  experiment  (June) 
on  days  one  (a),  two  (b),  three  (c),  and  four  (d). 
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Figure  A6.  Rhodamine  WT  distribution  in  Ranger’s  Slough  during  the  third  experiment  (June) 

on  days  five  (e),  six  (f),  and  seven  (g). 
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Figure  A7.  Rhodamine  WT  distribution  in  Bee  Creek  during  the  first  experiment  (February)  on 
days  one  (a),  two  (b),  three  (c),  four  (d),  and  five  (e). 
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Figure  A8.  Rhodamine  WT  distribution  in  Bee  Creek  during  the  second  experiment  (April)  on 
days  one  (a),  two  (b),  three  (c),  four  (d),  five  (e),  six  (f),  and  seven  (f). 
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Figure  A9.  Rhodamine  WT  distribution  in  Bee  Creek  during  the  third  experiment  (June)  on 
days  one  (a),  two  (b),  three  (c),  four  (d),  five  (e),  and  six  (f). 
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Figure  A10.  Rhodamine  WT  distribution  in  Thorp’s  Spring  during  the  first  experiment 
(February)  on  days  one  (a),  two  (b),  three  (c),  and  four  (d). 
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Figure  All.  Rhodamine  WT  distribution  in  Thorp’s  Spring  during  the  second  experiment  (April) 
on  days  one  (a),  two  (b),  three  (c),  four  (d),  and  five  (e). 
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Figure  A12.  Rhodamine  WT  distribution  in  Thorp’s  Spring  during  the  third  experiment  (June) 
on  days  one  (a),  two  (b),  three  (c),  four  (d),  five  (e),  and  six  (f). 
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Appendix  B:  Numerical  Modeling  Details 

This  appendix  presents  governing  equations,  variable  and  parameter 
definitions,  parameter  values,  forcing  parameters,  and  other  technical 
details  for  the  coupled  biological-physical  model  used  for  numerical 
simulations  of  P.  parvum  bloom  dynamics. 

Equations  for  the  biological  reactions  in  the  model  are: 


R z  PzZ  mzZ 

(B.la) 

Ri  =  PiNi~Zy,  i  —  p,c 

(B.lb) 

'  r 

Rr=  J2^iNi(iR4+ZYl  VR’1  VzQR'Z  ’  l=P’C 

i  i  Vi 

(B.lc) 

r  \ 

Rs=  vS>1  PzQs,z  ’  l=P’C 

'  i  {  Vi 

(B.ld) 

Rc=HcNcsc-kcC 

(B.le) 

Rr  =  {^'-^)NrSr-krP 

(B.lf) 

where 

Z  =  population  density  of  zooplankton 
Np  =  population  density  of  P.  parvum 
Nc  =  population  density  of  other  algae 
R  =  concentration  of  phosphate 
N  =  concentration  of  nitrate 

C  =  concentration  of  toxin  produced  by  cyanobacteria 
P  =  concentration  of  toxin  produced  by  P.  parvum 

Rz,  RP,  Rc,  Rr,  Rs,  Rc  and  Rp= respective  reaction  rates  for  the  above 
variables.  Table  Bi  lists  the  definitions  of  all  notation,  units,  and  parameter 
values  used  in  the  model.  The  equation  system  (B.i)  is  the  same  set  of 
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biological  reactions  and  interactions  used  in  a  previous  one-compartment, 
chemostat  model  for  Lake  Granbury  (Grover  et  al.  2012).  Further  explana¬ 
tion  of  the  assumptions  underlying  the  governing  equations  is  available  in 
this  publication. 


Table  Bl.  Notation,  units,  and  parameter  values 


Quantity 

Meaning 

Units 

Value 

State  Variables 

Zk 

Zooplankton  volume  density  in  compartment 
k;  k  =  1  for  main  lake,  k  =  2  for  cove 

pm3  liter1 

Variable 

Ni,k 

Population  density  for  algal  type  /  in 
compartment  k;  /  =  p  for  P.  parvum,  i  =  c  for 
cyanobacteria 

Cells  liter1 

Variable 

Rk 

Dissolved  phosphate  concentration  in 
compartment  k 

pmol  liter1 

Variable 

Sk 

Dissolved  nitrate  concentration  in 
compartment  k 

pmol  liter1 

Variable 

Ck 

Dissolved  cyanotoxin  concentration  in 
compartment  k 

pg  liter1 

Variable 

Pk 

Dissolved  P.  parvum  toxin  concentration  in 
compartment  k 

pg  liter1 

Variable 

Ecophysiological  Functions 

C  R  S 

Mz  ’  r^Z  ’  r'Z  ’  R’Z 

Zooplankton  population  growth  rate,  and 
potential  carbon-,  nitrogen-,  and  phosphorus- 
limited  growth  rates,  respectively 

d-1 

Variable 

fl 

Zooplankton  ingestion  rate  on  algal  type  / 

d-1 

Variable 

mz 

Zooplankton  mortality  rate 

d-1 

Variable 

Mi 

Population  growth  rate  of  algal  type  / 

d-1 

Variable 

//max,/' 

Nutrient-saturated  population  growth  rate  of 
algal  type  / 

d-1 

Variable 

/max 

Maximal  zooplankton  ingestion  rate 

d-1 

Variable 

Constant  Parameters 

ez 

Zooplankton  assimilation  efficiency  for 
volume  (carbon) 

None 

0.3 

Qr ,  Qs 

Phosphorus  and  nitrogen  quotas  per  unit 
zooplankton  volume,  respectively 

pmol  /  pm3 

0.061  x  10-9, 

1.2  x  10-9 

Qr.i,  Qs,/ 

Phosphorus  and  nitrogen  quotas  per  cell  for 
algal  type  / 

pmol  /  cell 

1.39  x  10-9,  32.0  x 

10  9  for  P.  parvum ; 
2.86  x  10-9,  32.5  x 
10-9  for 
cyanobacteria 

Vi 

Cell  volume  for  algal  type  i 

pm3  /  cell 

300  for  P.  parvum ; 

35  for  cyanobacteria 
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Quantity 

Meaning 

Units 

Value 

mz(  20) 

Zooplankton  mortality  rate  at  20  "C 

d-1 

0.14 

Qio 

Factor  by  which  zooplankton  mortality  and 
ingestion  rates  increase  for  a  10  °C  increase 
in  temperature 

None 

2.8 

Krj,  Ks,/ 

Half-saturation  constants  for  phosphorus- 
and  nitrogen-limited  growth  of  algal  type  /, 
respectively 

pmol  liter1 

0.009,  0.01  for  P. 
parvum ;  0.166, 

3.686  for 
cyanobacteria 

Kc 

Half-inhibition  constant  for  growth  of  P. 
parvum  in  response  to  cyanotoxin 
concentration 

pg  liter1 

0.1 

Exponent  controlling  nonlinearity  of  inhibition 
of  P.  parvum  growth  in  response  to  cyanotoxin 
concentration 

None 

1.0 

Wi 

Preference  for  zooplankton  ingestion  of  algal 
type  / 

None 

0.5  for  P.  parvum; 

0.49  for 
cyanobacteria 

Kz 

Half-saturation  constant  for  zooplankton 
ingestion 

pm3  liter1 

5.97  x  109 

KP 

Half-inhibition  constant  for  zooplankton 
ingestion  in  response  to  concentration  of  P. 
parvum  toxin 

pg  liter1 

0.5 

n 

Exponent  controlling  nonlinearity  of  inhibition 
of  zooplankton  ingestion  in  response  to 
concentration  of  P.  parvum  toxin 

none 

2 

fmax(20) 

Maximal  zooplankton  ingestion  rate  at  20  "C 

d-1 

2.4 

S’ 

Production  coefficients  for  cyanotoxin  and  P. 
parvum  toxin,  respectively 

pg  cell'1 

10.0  x  lO-s  for  P. 
parvum; 

12.5  x  10-9  for 
cyanobacteria 

kc,  kp 

First-order  degradation  coefficients  for 
cyanotoxin  and  P.  parvum  toxin,  respectively 

d-1 

0.5  for  both  toxins 

mT 

Optimal  growth  rate  of  P.  parvum 

d-1 

0.94 

i> 

Ratio  of  cove  volume  to  total  system  volume 
(cove  plus  main  lake) 

none 

Varies  for  cove 

¥ 

Ratio  of  flow  through  cove  to  total  system  flow 

none 

Set  equal  to  ^ 

E 

Exchange  rate  between  cove  and  main  lake 

d-1 

Varies  through 
range  bounding 
observed  values 

F 

Flushing  imposed  to  treat  blooms  of  P. 
parvum 

d-1 

Varies  according  to 
treatment  scenario 

ry  in  ry  in 

Zj  ,  z2 

Inflowing  zooplankton  density 

pm3  liter1 

Z[n  =  1  x  los; 

Z'"  =  0 
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Quantity 

Meaning 

Units 

Value 

Nin  Nin 

1  v  i,l>  1  v  i,  2 

Inflowing  algal  density 

Cells  liter1 

N‘]  =  0.1  x  106^ 

N‘"2  =  0 

n  in  n  in 

,  K2 

Inflowing  phosphate  concentration 

pmol  liter1 

PI"  =  R\"  =  forced 

o  in  o  in 

Oj  ,  o2 

Inflowing  nitrate  concentration 

pmol  liter1 

S[n  =  S2  =  forced 

s-iin  s~iin 

Lq  ?  L^2 

Inflowing  cyanotoxin  concentration 

pg  liter1 

C\n  =  C'"  =  0 

pin  pin 

r\  ’  r2 

Inflowing  P.  parvum  toxin  concentration 

pg  liter1 

P'n  =  If  =  0 

Forced  Parameters 

D 

Dilution  rate 

d-1 

Forced 

T 

Temperature 

°C 

Forced 

1 

Light 

pmol  rrr2  s1 

Forced 

G 

Salinity 

g  liter1 

Forced 

T)  in  o  in 

^  k 

Inflowing  phosphate  and  nitrate 
concentrations 

pmol  liter1 

Forced 

The  first  term  of  equation  B.ia  represents  growth  of  the  zooplankton 
population,  which  is  potentially  limited  by  the  rates  at  which  carbon, 
nitrogen,  or  phosphorus  are  acquired  by  consuming  algae.  The  per  capita 
growth  rate  pz  is  the  minimum  of  the  potential  growth  rates  limited  by 
these  elements: 


pz  =  min{pcz,pf,pFz} 


(B.2) 


where  the  potential  growth  rates  based  on  carbon  (C),  nitrogen  (N),  and 
phosphorus  (P)  consumption  rates  are  respectively  defined  as 

l4  =  ezY,fi  <B'3a> 


r  _  1  f&Rj 

v, 


,and  Vz  =  -prY 

Vs  i 


fiQs,i 


(B.3b,  c) 


In  these  expressions,  i  indexes  the  two  algal  types  represented  in  the 
model:  i  =  p  for  P.  parvum  and  i  =  c  for  cyanobacteria.  The  parameter  ezis 
the  efficiency  of  carbon  assimilation  by  zooplankton;  Qr  and  Qs  are  the 
phosphorus  and  nitrogen  quotas  per  unit  zooplankton  volume,  qRj  and  qs.i 
are  the  phosphorus  and  nitrogen  quotas  per  cell  for  algal  type  i,  and  Vi  is 
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the  cell  volume  of  algal  type  i.  The  function/-  is  the  ingestion  rate  of  algal 
type  i,  defined  in  detail  below  (Equation  B.7  and  following  text). 

The  second  term  of  equation  B.ia  represents  temperature-dependent 
mortality  of  zooplankton  with  per  capita  rate  mz : 


iT -20 

mz{T)  =  mz{20)Q[010 


(B.4) 


where 

T  =  temperature 
mz( 20)  =  mortality  rate  at  20°  C 

Q10  =  factor  by  which  zooplankton  mortality  increases  for  a  io°  C 
increase  in  temperature. 

The  first  term  of  Equation  B.ib  represents  growth  of  an  algal  population, 
either  P.  parvum  (i  =  p)  or  cyanobacteria  (i  =  c).  Algal  growth  rates  are 
potentially  limited  by  the  rates  at  which  nitrogen  or  phosphorus  are 
acquired  from  consuming  dissolved  nutrients.  The  per  capita  growth  rate  pi 
is  the  minimum  of  the  potential  growth  rates  limited  by  these  elements,  but 
it  is  also  temperature-  and  light-dependent  for  both  algal  types,  and 
salinity-dependent  for  P.  parvum.  These  complex  growth  dependencies  are 
represented  by  maximal  rates  of  algal  growth  that  depend  on  temperature, 
light,  and  salinity,  multiplied  by  the  minimum  of  Monod  saturation  terms 
that  depend  on  nutrient  concentrations.  For  P.  parvum,  an  additional  factor 
represents  inhibition  of  growth  by  the  dissolved  toxin  produced  by  cyano¬ 
bacteria.  Altogether,  the  per  capita  growth  rates  for  the  two  algal  types  are: 


PP=Pm^I>T^)mm 


R 


kr,p+R  ks,p+s 


Kl 


K(c+Ce 


(B.5a) 


Pc 


AW(7’r)min 


R 


kr,c  +  R’ks,c+s 


(B.5b) 


where 


//max,/  =  maximal  growth  rate  of  type  i 
I  =  irradiance 
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cr  =s  alinity 

KR,i  and  Ks,i  =  half-saturation  constants  for  growth  limitation  by  phosphorus 
and  nitrogen,  respectively 

Kc  =  half-inhibition  constant  for  the  effects  of  cyanotoxin  on  P. 
parvum 

t;  =  shape  parameter  controlling  the  shape  of  the  inhibition 
function. 


The  maximal  growth  rates  for  P.  parvum  and  cyanobacteria  are: 


p  ( a,I,T)  = 

•  max,p 

If  /  r-2o\| 

0,-3.531  +  0.02534-(ct  — 1.833)  +  7.468e  '  20  ^+0.000611(7-222) 

r  i T— 2°\i  (B.6a) 

-  0.06311  (a  - 1.833)2 -3.414e^4'  20  ^ -0.00000573(7 -222)' 

+0.1697(a- 1.833) 


/Wc U,D  =  rnax  0, 


1.683 1 


33.33  +  1  +  I2  /  300 


-0.01 


T- 5 
35 


(B.6b) 


Equation  B.6a  represents  a  unimodal  response  to  light,  temperature,  and 
salinity  for  P.  parvum  based  on  an  empirical  model  fitted  to  laboratory 
experiments  (Baker  et  al.  2009).  Equation  B.6b  represents  a  unimodal 
function  of  light  and  an  increasing  linear  function  of  temperature  for 
cyanobacteria  (Chapra  1997). 

The  second  term  of  Equation  B.ib  represents  mortality  for  both  algal  types 
due  to  zooplankton  grazing.  This  mortality  rate  depends  on  the  function/- 
for  the  consumption  rate  of  algal  type  i  by  zooplankton: 


fi(Np,Nc;T)  =  fmax(T) 


wyfN? 

Kp 

v  j 

Kp  +  p\ 

i,j  =  p,c  (B.7) 


The  consumption  rate  depends  on  algal  cell  densities  and  volumes 
according  to  a  sigmoid  functional  response  (Litchman  et  al.  2006),  with  a 
maximal  ingestion  rate  on  a  biovolume  basis  fmax  that  depends  on  tempera¬ 
ture.  The  parameter  Kz  is  a  half-saturation  constant  for  consumption  on  a 
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biovolume  basis,  and  the  parameters  uy  are  preferences  for  algal  type  i.  The 
consumption  of  algae  by  zooplankton  is  also  inhibited  by  toxin  produced  by 
P.  parvum,  through  a  factor  with  half-inhibition  constant  Kp  and  shape 
parameter  //.  The  temperature-dependence  of  zooplankton  consumption 
follows  the  function 


( T— 20 

fmJT)  =  fmax(20)Q^ 


(B.8) 


where 

/max(2o)  =  maximal  ingestion  rate  at  20°  C 

Qio  =  factor  by  which  consumption  rate  increases  for  a  io°  C 
increase  in  temperature. 

Both  algal  types  suffer  additional  losses  to  hydraulic  outflow,  which  are 
specified  by  the  transport  terms  in  the  physical  component  of  the  model. 
Other  loss  and  mortality  processes  are  neglected. 

The  first  terms  of  Equations  B.ic  and  B.id  for  nutrient  dynamics  represent 
consumption  of  dissolved  nutrients  coupled  to  algal  growth  through  the 
cell  quotas  for  nutrients.  The  second  terms  of  Equations  B.ic  and  B.id 
represent  recycling  of  nutrients  by  zooplankton  when  excess  nutrient  is 
acquired  through  ingestion  of  algae. 

The  first  term  of  Equation  B.ie  represents  toxin  production  by  cyano¬ 
bacteria,  proportional  to  their  productivity  through  a  coefficient  sc.  The 
second  term  of  Equation  B.ie  represents  first-order  decay  of  this  toxin  with 
rate  constant  kc.  The  first  term  of  Equation  B.if  represents  toxin  production 
by  P.  parvum,  which  is  proportional  to  the  reduction  of  growth  rate  below 
its  value  under  optimal  conditions,  with  coefficient  sp.  The  second  term  of 
Equation  B.if  represents  first-order  decay  of  this  toxin  with  rate  constant  kp. 

The  reaction  rates  represented  by  equation  system  B.i  are  placed  in  the 
physical  context  of  the  transport  equations  for  the  cove-main  lake  system. 
In  all  cases  the  reaction  rates  depend  only  on  the  values  of  variables  in  the 
respective  local  compartment,  cove  or  main  lake.  The  full  equation  system 
for  the  model  is  then: 
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dZ±  _  (1  —  </>)D  rym  D  +  E  v  ^D  +  £  +  F 
dt  l-i/;  1  1  - 1/;  1  l-i/; 


-Z2  +  Rz 


dZ2_<j>D^.in  E  v  <j)D  +  E  +  F  v 

—  z2  +  z1  Zj2-\- kz 

at  xp  xp  xp 


dN-‘  -  !1  '•’u’  V  .  _®±£w,.  +  *g±g±gy  +J^  i  =  p>c 


dt  1  -xp 


1-xp 


1  -xp 


^  =  tRNt  +  ENii_*D  +  E  +  FN 
at  xp  V  V 


dR,  (l-^)D  ;n  £>  +  £„  .  0D  +  E  +  F 


dt  1-xp 


-R™  - 


1-xp 


R i  + 


1-xp 


-R2  +  Rr 


dR2  _  <j)D  Djn  t  E  u  <pD  +  E  +  F 

—  K2  +  K± 

at  xp  xp  xp 


dS±  (1  —  </>)D  D  +  E  (j)D  +  E  +  F 


dt  1-xp 


-s;n- 


i-ip 


s .  + 


i-ip 


S2  +  Rs 


dS2  _  (j)D  cin  t  Ec  <j)D  +  E  +  F 

—  Oy 

dt  xp  xp  xp 


S2  +  Rs 


dC1  _  (1  —  </>)D  in  D  +  E  dD  +  E  +  F 

dt  1-xp  1  1-xp  1  1-xp 


-C2  +  Rc 


dC2  _  (f>D  nin  t  E  n  (pD  +  E  +  F 

7.  —  ^2  '  W 

dt  xp  xp  xp 


C2  +  Rc 


dP1  (1  -flDpin  D  +  Ep  xfrD  +  E  +  F  p 

1  a  _ . .  2  P 


dt  1-xp 


1-xp 


1-xp 


dP2  _  </)D  pm  E  p  tj)D  +  E  +  F 
dt  xp  2  xp  1  xp 


(B.9a) 

(B.9b) 

(B.9c) 

(B.9d) 

(B.9e) 

(B.9f) 

(B.9g) 

(B.9h) 

(B.9i) 

(B.9j) 

(B.9k) 

(B.9I) 


For  the  physical  model,  subscript  l  indicates  quantities  pertaining  to  the 
main  lake,  and  subscript  2  to  those  pertaining  to  the  cove.  The  parameter  (f> 
is  the  ratio  of  cove  volume  to  total  system  volume  (cove  plus  main  lake), 
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and  i// is  the  ratio  of  flow  through  the  cove  to  total  flow  through  the  system. 
The  parameter  D  is  the  dilution  rate  resulting  from  flow  through  the  whole 
system  (cove  plus  main  lake),  E  is  the  exchange  rate  between  cove  and 
main  lake,  and  F  is  the  additional  flushing  through  the  cove  imposed  by 
bloom  treatments.  The  parameters  with  superscript  “in”  are  inflowing 
concentrations  of  the  respective  constituents  into  the  cove  and  main  lake 
compartments.  Further  explanation  of  the  transport  parameters  of  the 
physical  model  and  the  values  adopted  in  this  study  are  provided  in 
Chapter  l  of  the  main  text  of  this  report. 

Six  of  the  parameters  in  the  model  are  time-variable:  dilution  rate  (D), 
salinity  (o),  irradiance  (7),  temperature  (T),  and  the  inflowing  concentra¬ 
tions  of  phosphate  and  nitrate  (7?f ,  S'" ).  These  parameters  were  changed 

daily  according  to  observational  data  as  compiled  in  Grover  et  al.  (2010). 
Daily  values  of  dilution  rate  were  obtained  by  estimating  Brazos  River 
discharge  as  the  mean  of  two  USGS  gauges:  08090800,  about  40  km 
upstream  of  the  headwaters  of  Lake  Granbury,  and  08091000,  about  60  km 
downstream  of  the  dam.  Discharge  was  then  divided  by  the  volume  of  the 
lake.  The  parameter  i//in  the  physical  model  partitions  this  flow  between 
cove  and  main  lake  compartments  and  was  assumed  to  equal  the  volume 
ratio  (f).  Other  data  were  obtained  from  approximately  monthly  sampling  of 
10  stations  at  Lake  Granbury  from  1  August  2006  to  31  August  2007,  which 
were  then  averaged  over  stations,  and  daily  values  were  calculated  by  linear 
interpolation.  These  daily  values  were  assumed  to  apply  equally  to  cove  and 
main  lake  compartments.  Nutrient  input  concentrations  were  assumed  to 
be  equal  to  measured  total  nutrient  concentrations  in  the  lake.  Average  daily 
irradiance  in  the  water  column  was  estimated  assuming  exponential 
attenuation  of  surface  irradiance,  with  an  attenuation  coefficient  estimated 
from  Secchi  depth,  for  an  average  depth  in  the  lake  of  5  m  (as  described  in 
detail  in  Grover  et  al.  2010). 

All  other  parameters  were  assigned  constant  values  (Table  Bi). 

Numerical  simulations  of  the  model  equation  system  (B.9)  were  computed 
with  an  adaptive  step-size,  fourth/fifth  order  Runge-Kutta  algorithm 
(Press  et  al.  1986).  This  algorithm  was  run  for  one  day  at  a  time,  using  the 
daily  values  of  the  six  time-variable  parameters.  The  model  was  coded  in 
Fortran  77  and  compiled  and  built  with  Absoft  Fortran,  version  1.0.4. 
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Appendix  C:  Dendrograms  and  Time-series 
Euclidean  Distances 


This  appendix  presents  dendrograms  and  time-series  Euclidean  distances 
for  the  cove-open  lake  station  comparisons  using  the  historical  data. 


Figure  Cl.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  A  and  the 
averaged  data  from  stations  2  and  3.  Here,  red  numbers  reflect  observations  from  the  cove 
and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired 
based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open 
lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake 
sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 
open  lake  sampling  during  April  2010. 
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Figure  C2.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic  (b) 
parameters  as  estimated  by  the  Euclidean  distance  for  comparison 
of  the  average  values  for  cove  A  and  stations  2  and  3. 
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Figure  C3.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  B  and  the 
averaged  data  from  stations  2  and  3.  Here,  red  numbers  reflect  observations  from  the  cove 
and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired 
based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open 
lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake 
sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 
open  lake  sampling  during  April  2010. 
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Figure  C4.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic  (b) 
parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  of  average  values  for  cove  B  and  stations  2  and  3. 
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Figure  C5.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  C  and  the 
averaged  data  from  stations  2  and  3.  Here,  red  numbers  reflect  observations  from  the  cove 
and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired 
based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open 
lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake 
sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 
open  lake  sampling  during  April  2010. 
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Figure  C6.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic  (b) 
parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  between  average  values  for  cove  C  and  stations  2  and  3. 
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Figure  C7.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  D  and  the 
averaged  data  from  stations  2  and  3.  Here,  red  numbers  reflect  observations  from  the  cove 
and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired 
based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open 
lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake 
sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 
open  lake  sampling  during  April  2010. 
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Figure  C8.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic  (b) 
parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  between  average  values  for  cove  D  and  stations  2  and  3. 
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Figure  C9.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  E  and  the 
averaged  data  from  stations  2  and  3.  Here,  red  numbers  reflect  observations  from  the  cove 
and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired 
based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open 
lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake 
sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 
open  lake  sampling  during  April  2010. 
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Figure  CIO.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic 
(b)  parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  between  average  values  for  cove  E  and  stations  2  and  3. 
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Figure  Cll.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  F  and  the 
averaged  data  from  stations  4  and  5.  Here,  red  numbers  reflect  observations  from  the  cove 
and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired 
based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open 
lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake 
sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 
open  lake  sampling  during  April  2010. 
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Figure  C12.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic 
(b)  parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  between  average  values  for  cove  F  and  stations  4  and  5. 
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Figure  C13.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  G  and 
the  averaged  data  from  stations  4  and  5.  Here,  red  numbers  reflect  observations  from  the 
cove  and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are 
paired  based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and 
open  lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open 
lake  sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 

open  lake  sampling  during  April  2010. 
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Figure  C14.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic  (b) 
parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  between  average  values  for  cove  G  and  stations  4  and  5. 
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Figure  C15.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  H  and 
the  averaged  data  from  stations  6  and  7.  Here,  red  numbers  reflect  observations  from  the 
cove  and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are 
paired  based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and 
open  lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open 
lake  sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 

open  lake  sampling  during  April  2010. 
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Figure  C16.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic 
(b)  parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  between  average  values  for  cove  H  and  stations  6  and  7. 
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Figure  C17.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  I  and  the 
averaged  data  from  stations  8  and  9.  Here,  red  numbers  reflect  observations  from  the  cove 
and  blue  numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired 
based  on  sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open 
lake  sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake 
sampling  during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and 
open  lake  sampling  during  April  2010. 
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Figure  C18.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic 
(b)  parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  between  average  values  for  cove  I  and  stations  8  and  9. 
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Figure  C19.  Dendrograms  comparing  biotic  (a)  and  abiotic  (b)  parameters  from  cove  J  and  the 
data  from  station  10.  Here,  red  numbers  reflect  observations  from  the  cove  and  blue 
numbers  reflect  observation  averages  from  the  open  lake.  The  numbers  are  paired  based  on 
sampling  month  and  year.  For  example,  [1,25]  represents  the  first  cove  and  open  lake 
sampling  during  December  2007;  [2,26]  represents  the  second  cove  and  open  lake  sampling 
during  January  2008;  and  so  forth,  with  [24,48]  representing  the  final  cove  and  open  lake 

sampling  during  April  2010. 
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Figure  C20.  Time-series  dissimilarity  between  biotic  (a)  and  abiotic 
(b)  parameters  as  estimated  by  the  Euclidean  distance  for  the 
comparison  between  cove  J  and  station  10. 
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